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Abstract 



The ANTARES Collaboration proposes to construct a large area water 

Cherenkov detector in the deep Mediterranean Sea, optimised for the detection of 
muons from high-energy astrophysical neutrinos. This paper presents the 
scientific motivation for building such a device, along with a review of the 
technical issues involved in its design and construction. 

The observation of high energy neutrinos will open a new window on the 
universe. The primary aim of the experiment is to use neutrinos as a tool to study 
particle acceleration mechanisms in energetic astrophysical objects such as active 
galactic nuclei and gamma-ray bursts, which may also shed light on the origin of 
ultra-high-energy cosmic rays. At somewhat lower energies, non-baryonic dark 

matter (WIMPs) may be detected through the neutrinos produced when 
gravitationally captured WIMPs annihilate in the cores of the Earth and the Sun, 
and neutrino oscillations can be measured by studying distortions in the energy 
spectrum of upward-going atmospheric neutrinos. 

The characteristics of the proposed site are an important consideration in detector 
design. The paper presents measurements of water transparency, counting rates 
from bioluminescence and potassium 40, bio-fouling of the optical modules 
housing the detectors photomultipliers, current speeds and site topography. 
These tests have shown that the proposed site provides a good-quality 
environment for the detector, and have also demonstrated the feasibility of the 

deployment technique. 

The present proposal concerns the construction and deployment of a detector 
with surface area 0.1 km^. The conceptual design for such a detector is 

discussed, and the physics performance evaluated for astrophysical sources and 
for neutrino oscillations. An overview of costs and schedules is presented. 

It is concluded that a 0.1 km^ detector is technically feasible at realistic cost, and 
offers an exciting and varied physics and astrophysics programme. Such a 

detector will also provide practical experience which will be invaluable in the 
design and operation of future detectors on the astrophysically desirable 1 km^ 

scale. 
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Chapter 1 
Introduction 



This document presents the scientific motivation for building a high energy 
neutrino undersea detector, with an effective area of 0.1 km^, along with a review 
of the technical issues involved in its design and construction. The planned 
apparatus consists of an array of photomultipliers, arranged in a lattice near the sea 
bed at a depth of 2400 m, to detect the Cherenkov light from muons produced by 
neutrino interactions in the seawater and rock beneath. Since 1996, the ANTARES 
collaboration has conducted an extensive R&D programme in close collaboration 
with experts in marine technology. This programme has addressed most of the 
critical technical aspects relevant to the construction and operation of a neutrino 
telescope, providing the information and experience necessary to be confident 
that the planned array is feasible at a reasonable cost. The neutrino detection 
capabilities of the design have been evaluated in detailed simulation studies and 
its potential science mission explored. 



1.1 Why neutrino astronomy? 

Most of our current knowledge of the Universe comes from the observation of 
photons. Photons have many advantages as cosmic information carriers: they 
are copiously produced, they are stable and electrically neutral, they are easy to 
detect over a wide energy range, and their spectrum carries detailed information 
about the chemical and physical properties of the source. Their disadvantage 
is that the hot, dense regions which form the central engines of stars, active 
galactic nuclei and other astrophysical energy sources are completely opaque 
to photons, and therefore we cannot investigate the properties of these regions 
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by direct observation, but only by indirect inference. For example, the photons 
we observe from the Sun come from its photosphere, far removed from the 
hydrogen-fusing core. Moreover, high energy photons interact with photons of 
the infrared radiation background and with the cosmic microwave background 
to create electron-positron pairs; this is the Greisen-Zatsepin-Kuz'min effect 
(GZK) [jT]]. This effect suppresses any possibility of surveying the sky over 
distances greater than 100 Mpc with high energy (>10 TeV) gamma rays. 

In order to observe the inner workings of the astrophysical objects and to 
obtain a description of the Universe over a larger range of energies, we need a 
probe which is electrically neutral, so that its trajectory will not be affected by 
magnetic fields, stable so that it will reach us from distant sources, and weakly 
interacting so that it will penetrate regions which are opaque to photons. The only 
candidate currently known to exist is the neutrino. 

Some astrophysical sources are known to emit neutrinos: hydrogen fusion 
produces electron neutrinos as by-products, and solar neutrino astronomy has a 
30 years long history; the conversion of iron nuclei to neutrons when a neutron 
star is formed in the heart of a supernova produces a burst of neutrinos (augmented 
by the thermal production of neutrino-antineutrino pairs), and one such burst was 
observed by Kamiokande and 1MB for Supernova 1987 A; cosmology predicts 
a low-energy relic neutrino background similar to the low-energy relic photons 
of the Cosmic Microwave Background, but these would have an effective 
temperature of around 1 .9 K and are very difficult to observe. 

Astrophysical sources of high-energy neutrinos have not been observed 
directly, but their existence can be inferred from the properties of cosmic rays. 
Primary cosmic rays are protons, with some admixture of heavier nuclei; the 
energy spectrum is a power law which extends to extremely high energies, values 
exceeding 10^° eV having been observed in recent years. Protons themselves have 
limited use as astrophysical information carriers because they are charged, and 
therefore subject to deflection by cosmic magnetic fields: only the very highest- 
energy cosmic rays are likely to retain any memory of the source direction. The 
exact source of the high-energy cosmic rays is thus unknown, although supernova 
remnants and active galactic nuclei have been proposed. Whatever the source, 
it is clear that accelerating protons to such high energies is likely to generate a 
large associated flux of photo-produced pions, which decay to yield gamma rays 
and neutrinos. These will remember the source direction, and so the existence 
of a general flux of very high energy cosmic -ray protons implies the existence of 
sources of high-energy neutrinos. 

Neutrino astronomy thus offers the possibility of observing sources which 



4 



correspond to the central engines of the most energetic astrophysical phenomena. 
As discussed below, it also provides long baselines for neutrino oscillation studies, 
and can explore useful regions of supersymmetric parameter space in the context 
of dark matter. The drawback, of course, is that the weak interactions of neutrinos 
imply that a very massive detector with extremely good background rejection is 
required to observe a measurable flux. 

1.2 The view from a neutrino telescope 

The ANTARES scientific programme is described in detail in chapter 2, but a brief 
overview is presented here. It is convenient to divide the programme into three 
broad subject areas: particle physics (neutrino oscillations), particle astrophysics 
(searches for neutralino dark matter) and astronomy. 

Within the minimal standard model, neutrinos are strictly massless, but 
the need to incorporate non-zero masses has long been anticipated. Neutrino 
astronomy was instrumental in establishing non-zero neutrino masses as an 
important topic: the lower- than-predicted flux of electron neutrinos from the 
Sun (the Solar Neutrino Problem) is now generally believed to be explained 
by neutrino oscillations either in vacuo or more probably within the Sun itself, 
enhanced by the high electron density in the solar core (the MSW effect). 

Neutrino oscillation solutions to the solar neutrino problem involve the 
conversion of electron neutrinos into some other flavour. The energies involved 
are not well matched to the detection capabilities of ANTARES. However, recent 
results from Super-Kamiokande appear to show a similar flux reduction occurring 
for muon neutrinos generated by cosmic-ray interactions at the top of the Earth's 
atmosphere. The effect is interpreted as evidence for the oscillation of into 
either ly^ or a so-called "sterile" neutrino, with a large mixing angle and a squared 
mass difference of ~ 10~^ — 10~^ eV^. This is a result of major importance which 
urgently requires confirmation by an experiment with independent systematics. 
As discussed below, the proposed ANTARES configuration should be capable of 
exploring the region of parameter space favoured by the Super-Kamiokande data. 

In recent years it has become generally accepted by astrophysicists that most 
of the matter in the universe is non-luminous "dark matter". The clearest evidence 
for this is the observed flatness of the rotation curves of disk galaxies, which imply 
a dynamical mass far in excess of that accounted for by the constituent stars and 
gas. Constraints from the observed abundances of light elements indicate that 
much of the dark matter in the cosmos must be non-baryonic. No presently known 
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particle has the required properties, but a good theoretical candidate is the stable 
neutral particle expected in most versions of supersymmetry theory. 

The detection and identification of a relic cosmological population of 
supersymmetric particles would be of immense importance to both cosmology 
and particle theory. Neutrino telescopes are not directly sensitive to a weakly 
interacting massive particle (WIMP). However, supersymmetric WIMPs will 
accumulate in the cores of the Sun and the Earth or in the centre of the 
Galaxy through gravitational capture. The resulting high space density leads to 
annihilation reactions, which will yield high-energy neutrinos through the decays 
of the gauge bosons and heavy particles produced. The proposed detector would 
be sensitive to these neutrinos over a useful range of WIMP masses. Compared 
to ongoing direct detection experiments, neutrino telescopes are generally more 
suitable for higher masses, although resonances in the Earth's capture cross- 
section enhance the signal strongly at certain lower masses, particularly around 
56 GeV. But even a confirmation of a prior direct detection would provide useful 
information about the couplings of the WIMP, and thus help to constrain the 
parameters of the theory. 

As discussed above, the two confirmed astrophysical neutrino sources are: the 
Sun and supemovae as exemplified by SN 1987 A. Both of these produce low 
energy neutrinos which could not be tracked by ANTARES, although a nearby 
supernova could be detected through a transient increase in the overall singles 
rate. However, there are a number of candidate astrophysical sources of high- 
energy neutrinos which would be detected by ANTARES. Since pions produced 
by high-energy protons are the likeliest source of high-energy neutrinos, candidate 
astrophysical neutrino sources are the proton accelerators, which might explain 
the High Energy Cosmic Rays spectrum. High energy 7-rays may be associated 
with high energy protons and their subsequent decays from 7r°, but they are also 
produced by synchrotron radiation of fast electrons in the presence of magnetic 
fields. The observation of neutrino sources would unambiguously discriminate 
between the two acceleration mechanisms. 

Candidate sources can be identified both within the Galaxy — accreting 
binaries containing neutron stars or black holes, supernovae and young supernova 
remnants — and elsewhere, most notably active galactic nuclei (AGN) and gamma- 
ray bursters. This represents a rich spectrum of possible sources, both steady and 
transient, covering a wide range of neutrino energies. 

It is worth noting that the history of astronomical observation suggests that 
the likeliest outcome of opening a new observational window is the discovery 
of a completely unexpected class of sources. Such discoveries are by definition 
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difficult to anticipate, but one possible pointer is that the very highest-energy 
cosmic rays, those above 10^° eV, remain difficult to explain in present models. 
This puzzle could be solved by the observation of their associated neutrinos. 

1.3 Present and future neutrino telescopes 

The proposed ANTARES detector is one of a number of present and proposed 
neutrino telescope projects. Of these, the Lake Baikal detector, consisting of eight 
strings supporting a total of 192 optical modules, was the first to demonstrate 
the feasibility of the technique, but is limited in depth and spatial extent by the 
nature of its site and so cannot be extended to a full-scale neutrino telescope. 
The AMANDA array at the South Pole uses ice rather than water as the detector 
medium, which gives it lower backgrounds but poorer directional accuracy. An 
array with an effective area of around 10^ m^, AMANDA B, has been deployed 
at a depth of 1500-2000 m and has published data on atmospheric muons; the 
installation of AMANDA-II, which will have an effective area several times larger, 
has started during the 1997-98 season. An eventual scale-up to a full km-scale 
array (ICECUBE) is proposed. Also at the proposal stage is NESTOR, with a 
deep site in the eastern Mediterranean off Pylos (Greece). 

The present ANTARES proposal is in many respects complementary to 
AMANDA. A water-based detector is more flexible and has better directional 
sensitivity, at the cost of higher background noise. Together, the two detectors 
cover the whole sky, with a substantial amount of overlap for cross-checking. 

The Mediterranean Sea represents an environment for a neutrino telescope that 
is quite different from those of the operating arrays at Lake Baikal (a freshwater 
lake which freezes over in winter) and the South Pole. Therefore, since its creation 
in 1996, ANTARES has followed a first phase R&D programme focused on three 
major milestones: 

• construction and deployment of test lines dedicated to measuring 
environmental parameters such as optical background, biofouling and water 
transparency; 

• development of prototype strings to acquire the necessary expertise to 
deploy and operate an undersea detector up to the kilometre scale; 

• development of software tools to explore the physics capabilities of the 
detector. 
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To date, the majority of the deployments have been performed in the 
Mediterranean Sea 30 km off the coast near Toulon (France), at a depth of 2400 m. 

This first phase has demonstrated that the deployment and physics operation 
of such a detector is feasible. The present proposal advances the ANTARES 
programme to a second phase, namely the construction of an array of 1000 optical 
modules to form a high-energy neutrino detector with effective area 0. 1 km^. Such 
an array would have the following achievable physics goals: 

• the study of the high energy neutrino flux (in the TeV-PeV range) with 
unprecedented angular resolution; 

• the measurement of atmospheric neutrino oscillations in the region of 
parameter space favoured by Super-Kamiokande; 

• a search for supersymmetric dark matter covering a region of model 
parameter space which is interesting for cosmology and particle physics. 

It would also provide practical experience and expertise which will be 
invaluable for the anticipated third phase during which a larger scale detector will 
be constructed, capable of conducting a search for astrophysical sources. 
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Chapter 2 

Scientific programme 



The science mission of a neutrino telescope such as ANTARES is vast, 
encompassing neutrino astronomy, neutrino oscillations and non-baryonic matter 
in the form of neutralinos, heavy metastable relic particles or topological defects. 
This chapter gives a brief account of the physics of the main candidate neutrino 
sources and the mechanisms by which they might produce a significant signal. 

2.1 Astrophysical sources 

The principal mechanism for generating high energy neutrinos is through decay 
cascades induced by high energy protons. Interactions of protons with matter or 

radiation produce mesons whose leptonic decay modes will yield neutrinos. A list 
of possible high energy astrophysical neutrino sources therefore involves a list of 
candidate astrophysical sites for proton acceleration. 

2.1.1 X-ray binaries 

This class of binary stars, which are among the brightest cosmic X-ray sources, 
consists of compact objects, such as neutron stars or black holes, which accrete 
matter from their normal companion stars. The accretion process leads to plasma 
waves in the strong magnetic field of the compact object, which bring protons to 
high energies by stochastic acceleration. Interactions of the accelerated particles 
with the accreting matter or with the companion star would then produce a 
neutrino flux comparable to that in high-energy particles with a spectral index 
close to 2. 
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2.1.2 Supernova remnants and cosmic rays 

Explosions of massive stars (supemovae) produce an expanding shell of material 
which is known from radio observations to accelerate high-energy particles. In 
some cases, the residue of the supernova is a neutron star which is detectable as a 
pulsar. Protons inside supernova shells can be accelerated by a first-order Fermi 
mechanism if (as seems likely) the shell is turbulent. If a pulsar is present there 
are additional acceleration mechanisms: in the magnetosphere of the pulsar, or 
at the front of the shock wave produced by the magneto-hydrodynamic wind in 
the shell [Q] . The interaction of these protons with the matter of the shell gives 
rise to neutrinos and photons (from charged and neutral pion decays respectively). 
An especially promising source is the ion doped wind model based on ultra high 
energy ion acceleration by the pulsar [^. The ions drift across the magnetic field 
in the surrounding supernova remnant filaments at a speed of order 10 times the 
hydrodynamic expansion velocity and interact with the thermal gas to produce 
neutrinos via n"^ production. The minimum flux predicted is substantially above 
background for the Crab remnant at 10^^ eV, but the theory of ion acceleration is 
too primitive to constrain the spectrum of injected ions and hence the production 
rate of energetic neutrinos at ANTARES energies. 

It is thought that supernova remnants are the principal galactic source of 
cosmic -ray protons, but, as yet, there is no direct confirmation of this hypothesis. 
Recent observations [Q| above 10^ eV by the EGRET detector have found 7— ray 
signals associated with at least 2 supernova remnants (IC 443 and 7 Cygni). 
However, electromagnetic radiation is an ambiguous signature, as it can come 
from either accelerated electrons or protons. An observation of neutrinos would 
provide a clear indication of proton acceleration with the direction identifying the 
source. 

2.1.3 Active galactic nuclei 

Active Galactic Nuclei (AGN) such as quasars are, averaged over time, the most 
powerful known objects in the Universe. Their observed total luminosities are 
in the range 10^^ - 10^^ W [|]]. In the generic model of AGN, these high 
luminosities arise from accretion of matter, at a rate of at least a few solar masses 
per year, onto a super-massive black hole, ranging from 10^ to 10^° solar masses. 
A minority of AGN also produces relativistic jets which transport synchrotron- 
emitting electrons, most easily observable in the radio, over distances up to 1 Mpc. 
Objects of this type whose jets are directed almost exactly towards us appear 
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as intense, compact and variable sources because their emission is amplified by 
Doppler boosting. Such objects are known as "blazars". 

An unexpected result, now well established by the EGRET satellite 
is that many blazars emit 7 rays with energies up to ~ 10^° eV. Even more 
surprisingly, four members of this class have been detected by atmospheric 
Cherenkov telescopes as highly variable emitters of 7 rays with energies 
exceeding 10^^ eV The mechanism by which such high energy particles are 
generated is a matter of active debate, as is the jet composition. 

Electrons are thought to be accelerated by the first-order Fermi process, and 
this is also efficient for protons. In generic AGN models, which do not necessarily 
involve jets, protons may be accelerated by shock waves associated with the 
accretion flow into the black hole. These protons may then interact with the matter 
of the accretion disk or with its ambient radiation field. Such models produce 
neutrinos with no associated high energy gammas, because the photons do not 
escape from the core of the AGN. 

In active galactic nuclei with jets, protons may also be accelerated in the 
inner regions of jets. These will then produce neutrinos by interaction with 
ambient radiation either emitted by the accretion disk or generated as synchrotron 
radiation within the jet. Similarly, Fermi acceleration in the hot-spots of powerful 
(Fanaroff-Riley Class II) radio sources may produce a population of high-energy 
protons [|]]. Such protons contribute to the high-energy 7-ray emission via a 
proton-induced cascade [Q, which will also generate neutrinos. Both neutrino 
and photon flux will be increased by Doppler beaming in the case of blazars. 

Recently limits have been set [|TD|] on the production of high energy neutrinos 
in AGN jets. These upper bounds, which are related to the measured flux of high 
energy cosmic rays, define fluxes of neutrinos significantly smaller than those of 
previous models. However, it remains possible that high-energy protons are found 
in some, but not all, AGN jets. 

Theoretical models of these sources are subject to large uncertainties, so that 
observations are of paramount importance. 

2.1.4 Gamma ray bursts 

Gamma Ray Bursts (GRBs) are the most spectacularly violent phenomena in the 
universe known to this date, and up until two years ago they were declared to be 
one of the outstanding mysteries of modern astrophysics. However, early 1997 
brought a major break though when the BeppoSAX satellite [ [TT[ ] located a burst 
precisely enough to permit the identification of its optical counterpart. Since 
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then, further advances have been made, with the discovery of GRB afterglows, 
measurements of their redshift and recognition of the magnitude of the energy 
release occurring in these second-long flashes. 

So far about a dozen GRB afterglows have been detected. The distances 
measured place the burst sources at cosmological distances with redshifts in 
the range z = 0.8 — 3.4 and indicate an energy release of 10"^^^^ J in 7- 
rays alone (assuming isotropic emission). The location of GRBs in their host 
galaxies somewhat correlates the population of GRB progenitors with the star 
formation rates, and supports models which involve a black hole formation 
through coalescence of a binary system of either a black hole-neutron star or a 
neutron star-neutron star. 

Detailed studies of afterglows from X-rays to optical and radio wavelengths 
provided crucial constraints on physical parameters for theoretical models. In 
the current standard model for 7-ray bursts and their afterglows, the fireball- 
plus -blastwave model [[T^, the initial event deposits a solar mass of energy into 
a region with a radius of about 100 km. The resulting fireball expands ultra- 
relativistically with Lorentz factor 7 > 300 into the surrounding medium. While 
the forward shock sweeps material and heats it, the reverse shock collides with 
the ejecta. Therefore, the afterglow is produced by synchrotron radiation when 
external shocks decelerate. However, the multi-peaked light curve of the gamma- 
ray burst itself is produced by the collisions of several internal shocks which are 
catching up each other with different Lorentz factors within the inner engine. 
There is clear evidence of these three distinct regions of the fireball in the light 
curve of GRB990123 [0. 

A hadronic component would naturally be expected in such extreme 
phenomena. Nearby GRBs could therefore be the long sought after sources of the 
highest energy cosmic rays. A possible example of such a relatively local GRB 
is GRB980425, whose location was observed to be coincident with the extremely 
bright supernova SN 1998bw. If this association is real, the GRB occurred at a 
distance of a few Mpc[|np and is a potential source of ultra-high-energy cosmic 
rays. In such a case, the cosmic ray observations set a model independent upper 
bound to the flux of high energy neutrinos produced by photo-meson interactions 
of the high energy protons with the radiation field of the source [[T0[]. This upper 
bound implies a muon flux from diffuse neutrinos of ~ 20/km^/year (in Att sr) in 
the energy range above 100 TeV. 

The rate predicted in reference [ [T0| ] is very low for ANTARES, which has 
a detector area of 0.1 km^ and an angular acceptance of 2ti sr. Nevertheless, a 
significantly higher event rate can be obtained by lowering the energy threshold 
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to around 100 GeV. The background is greatly reduced by requiring a spatial and 
temporal coincidence with an observed GRB, offering a unique opportunity for 
high energy neutrino detectors to observe neutrinos associated with individual 
bursts. 



2.1.5 Relic sources 

In recent years, the Fly's Eye atmospheric fluorescence detector and the AGASA 
air shower array have convincingly detected [ [T5| ] cosmic rays with energies 
exceeding ~ 5 x 10^^ eV — the Greisen-Zatsepin-Kuz'min (GZK) cutoff set 
by interactions on the 2.7 K black body cosmic microwave background. These 
ultra high-energy cosmic rays (UHECR) constitute a population distinct from 
those at lower energies (< 5 x 10^^ eV), in having a flatter spectrum. The 
depth in the atmosphere at which the shower reaches its maximum suggests a 
correlated change in the composition from iron nuclei to protons between 10^^ 
and 10^^ eV [[T6|]. The lack of any detectable anisotropy argues against a local 
origin in the Galactic disc where the presumed sources of low energy cosmic 
rays reside. However, there are no potential extra-galactic sources such as active 
galaxies near enough (within ~ 50 Mpc) to evade the GZK cutoff. Thus the origin 
of the UHECR is a major puzzle for standard physics and astrophysics. 

An exciting possibility is that UHECRs result from the decay of massive 
particles, rather than being accelerated up from low energies. The most 
popular models in this context are based on the annihilation or collapse of 
topological defects (TDs) such as cosmic strings or monopoles formed in the 



early universe OTTQ. When TDs are destroyed, their energy is released as massive 
gauge and Higgs bosons with masses of C(10^^) eV if such defects have formed 
at the GUT-symmetry breaking phase transition. The decays of these bosons can 
generate cascades of high energy nucleons, 7-rays and neutrinos. These models 



are constrained both by considerations of the cosmological evolution of TDs [ |T8| ] 



and observational bounds on the extra-galactic 7-ray background [ITPp. These 



require the mass of the decaying bosons to be less than ~ 10^^ eV, i.e. well below 
the GUT scale. Since only GUT scale TDs have independent motivation, e.g. to 
provide 'seeds' for the formation of large-scale structure, the above constraint thus 
disfavours TDs as the source of UHECRs. 

A more recent suggestion is that UHECRs arise from the decays of metastable 
relics with masses exceeding ~ 10^^ eV which constitute a fraction of the dark 
matter Such particles can be naturally produced with a cosmologically- 



interesting abundance during re-heating following inflation [Q. A lifetime 
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exceeding the age of the universe is natural if they have only gravitational 
interactions, e.g. if they are 'cryptons' — bound states from the hidden sector of 
string theory p3|, pT}]. This interpretation also naturally accounts for the required 
mass. A detailed study of the fragmentation of such heavy particles has been 
performed ^A\\ in order to calculate the expected spectra of nucleons, 7-rays and 
neutrinos. Such particles would, like all 'cold dark matter' particles, be strongly 
clustered in the Galactic halo, i.e. within ~ 100 kpc. Therefore, the extra-galactic 
contribution to the cosmic ray flux would be negligible in comparison and the 
observed flux fixes the ratio of the halo density to the lifetime. As an example, 
if such particles comprise all of the halo dark matter then the required lifetime is 
~ 10^° yr, with a proportionally shorter lifetime for a smaller contribution. 

Most of the energy in the cascade ends up as neutrinos and the predicted flux 
is then F(> E^,) ^ 10^{E^/1 TeVy^ km"^y"^sr"\ through normalisation to 
the observed UHECR flux. Furthermore, the neutrinos should be well correlated 
in both time and arrival direction with the UHECRs, given the relatively short 
propagation length in the halo. A small departure from isotropy of (9(20%) should 
also be observed, since our location is ~ 8 kpc from the galactic centre This 
anisotropy will be less than that for UHECRs, however, since there is no GZK 
cutoff to reduce the extra-galactic (isotropic) flux of neutrinos [^5]]. 

The essential point is that whatever process creates the UHECRs, it is 
exceedingly likely that there is a concomitant production of very high energy 
neutrinos. Measurement of the neutrino flux will, at the very least, provide 
important clues as to the origin of UHECRs and may even provide dramatic 
evidence for new physics. 



2.2 Neutrino oscillations 

The observation of atmospheric neutrinos with a neutrino telescope provides a 
means to study neutrino oscillations with a base-line length up to the order of the 
diameter of the Earth. The focus of investigation is the muon neutrino, but there is 
also a very interesting possible signature for extremely high energy (> 100 TeV) 
tau neutrinos. 



2.2.1 Atmospheric oscillations 

Atmospheric neutrinos are emitted in the decay of hadrons produced by the 
interactions of cosmic rays with atmospheric nuclei. The production of electron 
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neutrinos and of muon neutrinos is dominated by the processes ^ fi^ + I'^/v^ 
followed by fs^ ^ + V^/u^ + Vejve- In an infinite medium the ratio, r, of 
the flux of z/p and v^^ to the flux of z/g and V^. is expected to be two. Since the 
atmosphere is not an infinite medium, this ratio increases with increasing neutrino 
energy, because not all high energy muons can decay before they are absorbed 
by the ground. Furthermore, the magnetic field of the Earth has some influence 
on low energy charged particles and this modifies the energy spectra below a few 
GeV. 

These effects are taken into account for the calculations of the predicted 
neutrino fluxes. The overall normalisation uncertainty is estimated to be about 
20%, which is due to systematic theoretical uncertainties in the energy spectra of 
the primary cosmic rays and to the uncertainties in their composition. Generally, 
experimental results are reported as i? = tdata/'^mc in order to cancel common 
systematic uncertainties thereby reducing the overall uncertainty in R to about 
5%. 



Measurements published by underground experiments |2^, g9p show 
evidence of a deficit in the number of muon neutrinos with respect to electron 
neutrinos. No anomaly was observed by the Frejus [ |30| , |3T| ] and NUSEX [ ]32| ] 
experiments. 

Neutrino oscillations have been proposed as an explanation for the low value 
of the ratio R. With the hypothesis of two-neutrino mixing, the oscillation 
probability is: 



P = sin^ 29 sin^ (^1.27^Am'^^ 

where 9 is the mixing angle, L is the distance travelled by the neutrino (in km), 
E is the neutrino energy (in GeV) and Am^ is the difference of the square of the 
masses (in eV^). As the neutrinos are produced in the atmosphere, the distance 
L ranges between 15 km, for vertically downward-going neutrinos, and almost 
13 000km, for vertically upward-going neutrinos. 

A recent analysis [ ]2^ ] reported by the Super-Kamiokande collaboration 
outlines evidence for ^ u^, where Ux may be z/^ or a "sterile" neutrino, 
with sm^29 > 0.82 and 10"^ < Am^ < 8 x IQ-^eV^ at 90% confidence 
level. The most probable solution is Am? = 3.5 x 10~^eV^ and sin^ 26* = 1.0 
(maximum mixing). Given these values, the survival probability 1 — P vanishes 
for L/E = (2n+l) x 353 km GeV^^, where n is a non-negative integer. Figure [2?T| 
shows the variation of the survival probability as a function of L/E. 



15 




10 10 
L/E (km/GeV) 



10 



Figure 2.1: Variation of the neutrino survival probability as a function of L/E for 
Am^ = 3.5 X lO^^eV^ and maximal mixing. 



Several points remain to be clarified. First of all, there are three different 
regions of Am? which have been explored with solar neutrinos, atmospheric 
neutrinos and short baseline beam neutrinos respectively, all of which indicate 
evidence of non-zero Am^. This cannot be supported in a three-flavour neutrino 
scheme. Moreover, the present result found in the Super-Kamiokande experiment 
is only marginally compatible with the result of Kamiokande [ |2^ ] (Am? = 1.8 x 
10^^) and is in disagreement with the analysis of stopping upward-going muons 
performed by the 1MB collaboration [ pB] ] which excluded the region presently 
favoured by Super-Kamiokande on the (sin^ 29, Am^) diagram. The results 
reported by the Frejus [ PD[ ] experiment exclude the region Am? > 4 x 10^^ eV^ 
at 90% confidence level. 

An independent measurement of the region 60 < L/E < 1250km/GeV could 
resolve the uncertainties in the interpretation of these different experiments. This 
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region contains the principal oscillation (first dip) for the Am^ values in question: 
1 X 10"^ < Am^ <2 X 10^2 eV^. ANTARES is perfectly suited to this task. 
The peak sensitivity of the ANTARES detector is very near the most probable 
value reported by Super-Kamiokande. Accurate measurement of the position of 
the principal oscillation would provide a precise measurement of Am^. 

Determination of the neutrino survival probability shown in figure [2?T| requires 
a measurement of the energy of the incident muon neutrino. For the isotropic 
vq and T7q charged-current interactions, half of the neutrino energy goes to the 
hadron shower. For the T7q and I'q interactions, a much larger fraction of the 
energy goes to the muon, but these interactions are three times less frequent 
than the isotropic interactions. The energy of the hadron showers is difficult 
to estimate accurately, so the measurement of the oscillation parameters in 
ANTARES depends principally on the measurement of the muon momentum. The 
energy of the muon is determined by its range and the precision depends on the 
vertical spacing of the PMTs. The scale is set by the energy loss for vertical muons 
between groups of PMTs: 2 GeV for a vertical spacing of 8 m, and 4 GeV for a 
spacing of 16 m. Reconstruction inefficiencies degrade the energy resolution in 
an energy-dependent way, depending on the reconstruction algorithms employed. 
The current status of the reconstruction effort is described in chapter ^. 



2.2.2 Tau neutrinos 

Although the contribution of tau neutrinos to the atmospheric neutrino flux is 
negligible, their interactions must be taken into account when studying neutrino 
oscillations. If oscillations of the type z/^i Ur occur, the charged current 
interactions of u^. will produce charged r leptons that can contribute to the signal 
observed. Most of the hadronic and electronic r decays will escape detection, but 
the muonic decays fi^u^Ur can be seen in the ANTARES detector, and 

these events might be mis-identified as interactions. 

The muonic branching ratio of the r is 17%. Furthermore, the charged-current 
interactions are considerably suppressed by the limited phase space due to the 
mass of the r (1.78 GeV). For very large Am^ (1 eV^), the number of would 
be equal to the number of z/^ over the entire atmospheric flux (< 500 GeV), 
but the number of r^(r+) produced with energies above 10 GeV would be only 
53% (64%) of the number of because of the limited phase space, and the 

contamination of the muon sample due to tau decays would be about 9% (11%). 
For smaller values of Am?, the contamination in the region of the main oscillation 
dip could be larger. 
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In P^ ] an attractive possibility for detecting very high energy neutrinos is 
proposed. The Earth is nearly transparent to low-energy neutrinos, but opaque 
to neutrinos above 100 TeV. Nonetheless, tau neutrinos well above 100 TeV can 
produce a signal in ANTARES because, unlike the e"^ and fi"^ produced in z/g and 

interactions, the t"^ produced in z/^ interactions decay before they are absorbed, 
producing z/^ of lower energy which continue along the original flight path, but 
with decreasing interaction probability. Once the energy has been degraded to 
about 100 TeV, the u^. can penetrate the Earth and produce an accumulation of 
very-high-energy events in the detector. Such an accumulation would be a signal 
for tau neutrinos. Moreover, the flux of Ur from a given source would be constant 
during the Earth's rotation, whereas the flux of would vary with the sidereal 
day because of the change in elevation seen from ANTARES. The variation of the 

flux could not be observed from AMANDA, because it is located at the South 
Pole. A comparison of signals from the same sources observed at the different 
elevations corresponding to AMANDA and ANTARES could lead to very exciting 
results. 



2.3 Indirect detection of neutralinos 
2.3.1 Neutralinos as dark matter candidates 

There is extensive astrophysical evidence [ P5| ] that most of the matter in the 
universe is non-luminous. The matter content of the universe is normally 
described in terms of the density parameter ^7 = p/pc, where pc = SHq/StiG 
is the critical density and Hq, often expressed in the dimensionless form h = 
Hq/IOO km s^^ Mpc^^, is the expansion rate. Estimates of galactic halo masses 
from rotation curves and of the galactic number density give a value of > 0.1, 
whereas the luminous matter in galaxies corresponds to ^7 < 0.01. Studies of 
the dynamics of clusters and superclusters of galaxies increase the required value 
of ri to around 0.2-0.3, and the popular inflationary paradigm generally requires 
= 1 (although this may include a contribution from a non-zero cosmological 
constant). 

The density of baryonic matter is constrained by the abundances of the light 
elements to a range 0.008 < fibaryons^^ < 0.024. With current measurements 
of h tending to lie in the range 0.6-0.7, this indicates that much of the dark 
matter required on the galactic scale and beyond must be non-baryonic. A prime 
candidate for non-baryonic dark matter is the Lightest Supersymmetric Particle 
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(LSP). 

Supersymmetry is a spontaneously broken symmetry between bosons and 
fermions, postulated as a natural mechanism for avoiding GUT-scale radiative 



corrections to the Higgs mass. The minimal supersymmetric model (MSSM) pq ] 
contains boson partners for every fermion, fermionic partners for all known 
bosons, and two Higgs doublets. In order to match the experimental limit on 
proton decay, it is generally assumed that a multiplicative quantum number 
called R-parity is conserved: R = +1 for ordinary particles and —1 for their 
supersymmetric partners. The natural consequence of this is that the lightest 
supersymmetric particle is stable, having no R-conserving decay mode. 

The MSSM has numerous free parameters (105 from the most general soft- 
breaking terms in addition to the 19 parameters of the Standard Model, reducing 
to seven if we make a number of standard — but not necessarily correct — 
assumptions about masses and mixing angles), and the identity of the LSP is not 
unambiguous. However, experimental limits rule out a charged LSP over a broad 
mass range p7|], and it is therefore generally assumed that the LSP is the lightest 
of the four neutralinos, the mass eigenstates corresponding to the superpartners of 
the photon, the Z and the two neutral Higgs bosons. The mixing which generates 
the mass eigenstates depends on the choice of SUSY parameters: the LSP, x, 
can be anything from a nearly pure B-'mo to a nearly pure higgsino P^]. The 
mass of the LSP is constrained from below by non-detection in LEP2, and from 
above by the requirement that supersymmetry fulfil its role of maintaining the 
mass hierarchy between the GUT scale and the electroweak scale. 

Neutralinos produced in the early universe will fall out of equilibrium when 
the annihilation rate, (cr^f )n^, falls below the expansion rate H. This condition 
leads to the estimate that 

3 X 10^27 ^jj^3 



"x — 



where {aAv) is the thermally averaged annihilation cross-section times the relative 
velocity. (Note that neutralinos are Majorana particles, x = X-) If the masses of 
SUSY particles are close to the electroweak scale, the annihilation cross-section 
is of order (100 GeV)^ ~ 10"^^ cm^ s"^, indicating that ~ 1 is a realistic 
possibility. This heuristic argument is confirmed by detailed calculations, which 
indicate that a neutralino LSP has a cosmologically significant relic abundance 
over large regions of the MSSM parameter space [F^ 
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2.3.2 Neutralino detection 



If neutralinos make up a significant fraction of the Galactic dark halo, they will 
accumulate in the core of bodies such as the Earth or the Sun [ |38t |39| , ^ and in the 
Centre of our Galaxy pTl). A neutralino passing through such a body has a small 
but non-zero probability of scattering off a nucleus therein, so that its velocity after 
scattering is less than the escape velocity. Once this happens, repeated passages 
through the body will generate additional scatters and the neutralino will sink 
relatively rapidly to the centre. Equilibrium will be reached when the gain of 
neutralinos from capture is balanced by the loss from annihilation, C = CaN"^ 
where C is the capture rate, N is the number of captured neutralinos. The quantity 
Ca depends on the WIMP annihilation cross-section and the WIMP distribution. 

The capture rate C depends on the local halo mass density p^, the velocity 
dispersion of neutralinos in the halo v, and the elastic scattering cross-section, 
which depends on the mass of the neutralino and on the effective volume and 
chemical composition of the Sun or Earth. Representative values of and v are 
0.3 GeV cm^^ and 270 km s^^ respectively. Jungman et al. [ |38| ] quote capture 
rates of 

Cq = 2.4 X 1037 s-i ^/0(mJ/p(GeV')2, 

° 0.3 GeV cm-^'^^^ '''■"^^ ' ' 

= 2.4 X 10'^ s-^ ^ s/®(^x)/p(GeV')^ 

® 0.3 GeV cm-^'^®^ '"•'^^ ' ' 

where fp is the neutralino-nucleon scalar coupling and the dependence on the 
neutralino mass is given by /© and /® for Sun and Earth respectively. For 
the Sun, is a smoothly varying function of mass, decreasing from 0(5) at 
= lOGeV/c^ to C(0.05) at 1 TeV/c^; for the Earth, resonances occur 
whenever the neutralino mass equals the mass of the target nucleus, so there is 
a good deal of structure in below = 100 GeV/c^, with a main resonance at 
56 GeV (56Fe). 

The annihilation rate can be calculated [^, ^ from the thermally-averaged 



cross-section, (ctav), in the limit of zero relative velocity; at equilibrium it is 
simply half the capture rate. The time taken to reach equilibrium, r = (CC^)"^^^, 
is longer for the Earth than for the Sun: over a wide range of MSSM parameter 
space, the Earth's neutralino population has not yet reached equilibrium, and the 
observed annihilation rate will be suppressed by some model-dependent amount 
compared to the equilibrium rate. 
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2.3.3 Neutrinos from neutralino annihilation 



Neutralinos annihilate into a fermion-antifermion pair or into various two-body 
combinations of W, Z and Higgs bosons Direct decay into neutrinos is 

zero in the non-relativistic limit, but decays into c, b and t quarks, r leptons, 
Z, W and Higgs can all produce a significant flux of high-energy neutrinos 
(light quarks and muon pairs do not contribute, as they are stopped before they 
decay). The typical neutrino energy produced is thus around one-half to one-third 
of the neutralino mass, with a broad spectrum whose detailed features depend 
on the branching ratios into the different channels (which in turn depend on 
the neutralino composition — gaugino vs higgsino — as well as its mass) and are 
modified, especially in the Sun, by hadronisation and stopping of c and b quarks, 
and stopping, absorption and possibly oscillation of neutrinos ^ 

Since both the neutrino-nucleon cross-section and the range of the produced 
muon scale with the neutrino energy, the resulting muon rate is approximately 
proportional to dN^/dE^, indicating that this method of neutralino detection 
is most likely to be competitive for higher mass neutralinos. 

The expected muon flux depends on the neutralino mass and the assumed 
MSSM parameters. The results of representative calculations [ |39| ] are shown in 
figure In these models, the flux from the Sun exceeds that from the Earth at 
high neutralino masses (in others, e.g. [pO|], the solar flux is larger throughout). 



The exposure required to detect a signal (defined in [ |39| ] as a 4a effect above 
background with at least 4 events) is of the order of 10^ to 10^ m^-yr. 

The Sun is a point-like source of neutrinos from neutralino annihilation, but 
the Earth is not, especially for lower mass neutralinos. Despite the smearing 
produced by measuring the muon flux rather than the neutrinos, a detector with 
good angular resolution would be able to use the observed angular distribution to 
constrain the neutralino mass, as shown in figure ^3] from [pO|]. 

Indirect detection explores a region of parameter space which is somewhat 
different from that investigated by direct techniques, although there is also 
considerable overlap. Indirect detection is particularly useful when the axial 
coupling of the neutralino dominates pS]], as axially-coupled neutralinos are 
captured in the Sun due to their interaction with hydrogen. Detection of 
neutralinos by both direct and indirect techniques would provide information on 
neutralino couplings and MSSM parameters, as would indirect detection of signals 
from both the Sun and the Earth. 
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Figure 2.2: Exposure required to detect a neutralino signal. The figure is taken from 
Bottino et al. [39] and the 3 curves correspond to different choices of supersymmetry 
parameters. 



2.4 Other exotic phenomena 
2.4.1 GUT monopoles 

Grand Unified Theories (GUT) of the electroweak and strong interactions predict 
the existence of massive magnetic monopoles. It is assumed that they were 
produced shortly after the Big Bang and would have cooled down to very low 
velocities (/5 < 10^^). They may be at the origin of baryon number violating 
processes and have been searched for in proton decay experiments. 

The best limits to date are provided by the MACRO and BAKSAN [|5p 
experiments. In the region 10"'^ < /9 < 10"'^, they have achieved sensitivity below 
the Parker bound p6|] . Baikal has obtained limits which could be comparable with 



a limited number of photomultipliers (36 and 192) [gj. A 1000 PMT detector 
could bring appreciable improvement on the results. ANTARES is studying a 
special trigger similar to that of Baikal. The slow moving monopole produces 
sequential Cherenkov flashes along its track via the proton decay products. The 
trigger is based on counting rate excesses during short time periods of 100 
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Figure 2.3: Mass of the neutralino versus expected angle for different supersymmetry 



parameters. The figure is taken from Edsjo [40]. 



microseconds to 1 ms. 



2.4.2 Relativistic monopoles 



According to [p^, the Ultra High Energy Cosmic Ray Events above 10^° eV may 
be relativistic monopoles. In this case, the basic mechanism for light generation 
is Cherenkov radiation. The intensity of the light due to relativistic monopoles 
is giant, similar to that of PeV muons, and they can be obtained with the muon 
trigger. Preliminary estimates show that they can be distinguished from high 
energy muons, down to monopole velocities /3 = 0.1. 
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2.4.3 SUSY Q-balls 



Super symmetry predicts the existence of non-topological solitons, often called Q- 
balls. It has been realized that solitons with a large number of baryons are stable 
and can be produced copiously in the early Universe, which means that they are 
an interesting dark matter candidate [|49|]. The signature and rate of occurrence of 
a Q-ball through a detector depend on the parameters of the theory. In the case of 
large baryon numbers of the order of 10^'', emissions of the order of 10 GeV/mm 
are possible, a spectacular signature, similar to those expected from monopole 
searches [^. 

2.4.4 New phenomena 

Whenever a new window is opened on the Universe, unexpected phenomena are 
observed. This has been verified in numerous cases throughout history, from 
the observation of the moons of Jupiter by the first Galilean telescopes to the 
discovery of pulsars with radio-based astronomy. Neutrino astronomy provides an 
exciting new probe of the Universe to the highest possible red-shifts in an energy 
window so far not observable by other techniques. There is plenty of scope for 
surprises. 
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Chapter 3 

Detection principles 



Since the Earth acts as a shield against all particles except neutrinos, a neutrino 
telescope uses the detection of upward-going muons as a signature of muon 
neutrino interactions in the matter below the detector. The muon detection 
medium may be a natural body of water or ice through which the muon emits 
Cherenkov light. Its detection allows the determination of the muon trajectory. 
This detection technique requires discriminating upward going muons against the 
much higher flux of downward atmospheric muons (figure pTTt ). To simplify the 
discrimination, the detector should be installed in a deep site where a layer of 
water or ice would shield it. 

In order to correlate the measured muon spectrum with the original neutrino 
spectrum, it is necessary to understand the dynamics of neutrino interactions, the 
opacity of the Earth, the energy loss of muons and the resolution of the detector 
over a wide range of angles and energies. 



3.1 Neutrino interactions 

The inclusive deep inelastic charged current cross-section for ui + N ^ + X 
(where the lepton mass is neglected) is given by [ J5T1 ] 



xq{x,Q^)+x{l-yYq{x,Q^) 



dxdy TT (g2 + M^) 

where is the Fermi constant, m^v and Mw are the nucleon and 14^ -boson 
masses, and is the square of the momentum transfer between the neutrino and 
muon. The Bjorken variables x and y are x = Q"^ /2miqv and y = {Ey — Ei)/Ey, 
where u = E^ — Ei h the lepton energy loss in the laboratory frame. 
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Figure 3.1: Zenith angular distribution of the muon flux above 1 TeV from atmospheric 
muons and atmospheric neutrino induced muons at 2300 m water equivalent depth. 



3.1.1 Interactions at low energies 

At energies such that E^^ <^ M^/2m7v ~ 5 TeV, can be neglected in the 
W propagator. In this case the average deep-inelastic uN cross-sections grow 
linearly with the neutrino energy (see figure 

cTDisii^N) ~ 0.67 X IQ-'^^E^GeV] cm^ 
(TDisli^N) ^ 0.34 X lO-^s^^GeV] cm^ 

At even lower energies {Ei, < 100 GeV), quasi-elastic and resonant 
contributions to the charged-current interactions have to be considered. In this 
case Q"^ has to be small enough to allow a coherent interaction with the complete 
target nucleon, so these cross-sections are essentially constant with energy. A fit 
to data above 10 GeV gives the following sums for the quasi-elastic and resonant 
production processes [Q: 

<yQE{i^f.N) + ares{i^^N) = 1.50 X lO-^^ cm^ 
(Tqe{^i.N) + aresii^^N) = 1.58 X 10-38 cm2 

Between 10 GeV and 500 GeV the vertical z/^ flux is approximately proportional 
to l/Ef,. Integrating the cross-sections given above over this flux for .Ej, > 10 GeV 
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Figure 3.2: Average charged-current cross-section for z/-N interactions for different 
sets of paiton distribution functions. The data point coiTcsponds to the average of the 
measurements by HI and ZEUS at HERA (taken from [^]). 



results in a quasi-elastic and resonant contribution of 1 1% for z/^ and 21% for t'^ 
interactions. 



3.1.2 Interactions at high energies 

At energies such that Ey ^ M^/2mj\f ^ 5 TeV, quasi-elastic and resonant 
contributions are completely negligible. On the other hand, the W propagator 
limits the growth of to {Q"^) ~ and so the cross-section is dominated by 
the behaviour of distribution functions at small x. 

The HI and ZEUS collaborations at HERA measured the proton structure 
function F2{x, Q"^) via charged current e — p scattering, for in the range from 
1.5 to 5000 GeV^ with x down to 3 x IQ-^ at = 1.5 GeV^ and x down to 
2 X 10-2 at g2 = 5000 GeV^ [|5|]. 



These measurements can be translated into a neutrino-nucleon interaction 
cross-section at ~ 50 TeV and can also be used as a guide to extrapolate 
the parton densities beyond the measured ranges in x and to those required for 



higher neutrino energies. Figure shows the behaviour of the average vN cross- 
section for different sets of parton distribution functions. At very high energy, 
the cross-section calculated with the new parton distribution functions (CTEQ3- 
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DIS [|5^) is more than a factor of 2 larger than previous estimates. 



3.1.3 Different types of v interactions in ANTARES 

Charged-current z/g interactions give rise to electromagnetic and hadronic showers 
with longitudinal dimensions of no more than a few metres, because the radiation 
length and the nuclear interaction length of water are below 1 m. On the scale of 
ANTARES, these are nearly point-like events. At energies above 100 GeV, the 
energy resolution of these events is expected to be better than for muonic events 
because they leave all of their energy inside the detector volume. On the other 
hand, their angular resolution will be poor compared to muonic events, due to the 
point-like character of the showers. 

Charged-current interactions will be contaminated by neutral-current 
interactions of both z/g and z/^j (and z/^, if present). The number of neutral- 
current interactions is about 1/3 of the number of charged-current interactions. 
The neutrino type is not identified in the neutral-current interactions, the energy 
resolution is poor due to the missing final-state neutrino, and the angular 
resolution is poor due to the point-like character. 

Charged-current interactions produce fi^ leptons as well as a point-like 
hadronic shower. The z/^ energy can be estimated from the measured /i^ energy. 
In u^d jjL^u interactions, the average /i^ energy is 1/2 of the energy; 
in u^u fi^d interactions, the average jj,^ energy is 3/4 of the z/^ energy. 
The js^ energy can be determined from the range for < 100 GeV, or from 
dE/dx for E^ > 1 TeV (see below). For z/^ interactions inside the detector, 
additional information on the z/^ energy is available from the hadronic shower. 
The ANTARES detector is designed for the detection of these charged-current z/^ 
interactions. 

Charged-current Ur interactions produce leptons with electronic, muonic 
and hadronic decay modes. The Ur interaction vertex and the decay vertex 
cannot be separated for energies below ~100 TeV. The electronic and hadronic 
modes will look like z/g charged-current or neutral-current interactions. The 
muonic decays yU^z/^jZ/,-, with branching ratio 17%, will be visible in 

ANTARES, but they cannot be distinguished from z/^ interactions. 
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3.2 Cherenkov light emission 



Charged particles emit light under a characteristic angle when passing through a 
medium if their velocity exceeds the speed of light in the medium. The Cherenkov 
angle 9 is related to the particle velocity /3 and the refractive index of the medium 
n: 

cose^l/n(3 (3.1) 

In the energy range interesting for ANTARES (E > 10 GeV), particles will 
generally be ultra-relativistic with P — 1. The refractive index of sea water is 
n = 1.35 for a wavelength of 450 nm therefore the Cherenkov light is emitted 
under 42° for this wavelength. This easy geometrical pattern of light emission 
allows a precise reconstruction of tracks from the measurement of only few hits 
at different space points. 

The number of photons produced along a flight path dx in a wave length bin 
dX for a particle carrying unit charge is 

= 27rQ;sin2^/A^ (3.2) 

dXdx 

At wavelengths of 400-500 nm the efficiency of the photomultipliers as well as the 
transparency of the water are maximal. Within 1 cm flight path 100 photons are 
emitted in this wavelength bin. Between 285-400 nm twice as many photons are 
emitted, however they contribute less to the detected signal. At a perpendicular 
distance of 40 m from a charged track the density of photons between 400- 
500 nm is still 1 per 340 cm^, neglecting absorption and scattering effects. The 
effective area of the photomultipliers being considered is in the same range (300- 
500 cm^). This gives an indication of the active detector volume around each 
photomultiplier. 

For P — 1 the Cherenkov light yield is independent of the energy of the 
charged particle. This means the light output of a single particle does not 
allow its energy to be measured. However when hadronic or electromagnetic 
showers are produced (which might occur at the neutrino vertex as well as for 
radiative processes along a muon track) the total light yield of the shower will 
be proportional to the total track length in the shower and therefore to its initial 
energy. This allows some calorimetric measurements if the neutrino vertex is 
inside the active detector volume or for muon tracks above 1 TeV where radiative 
processes dominate its energy loss. 
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3.3 Light propagation in sea water 



The processes of absorption and scattering characterise the transmission of light 
in water. They are parametrised by the absorption length A^, the scattering length 
\s, and the scattering function P(9) which describes the angular distribution of the 



scattering [ p5| ] . The relevant window of wavelengths for a sea water Cherenkov 
detector is centred on blue light. Deep sea water transparency is maximal in the 
blue, with typical values of 60 m for \a and A^, and a scattering function peaked 
in the forward direction with an average value for the cosine of the scattering 
angle (cos(^)) 



0.9, as shown in figure |3J| taken from [p5||. Seasonal variations 
are expected to affect these values, especially the scattering parameters which 
are governed by the amount of suspended particulate matter. In section |4.2.3| we 



present results of in situ measurements of optical properties at the ANTARES site. 
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Figure 3.3: The angular distribution of scattering in the deep Mediterranean Sea [pSQ. 



3.4 Detector response 

The physical processes involved in neutrino and muon interactions place limits on 
the angular and energy resolution possible with a neutrino telescope. These limits 
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must be taken into consideration when optimising the detector design. 



3.4.1 Angular response for interactions 

The angular response of the detector with respect to the incoming neutrino 
direction is crucial for the identification of point sources of neutrinos. Three 
factors determine this response: the angle between the neutrino and the muon 
in the neutrino interaction, the deviation of the muon direction due to multiple 
scattering and the angular resolution of the detector with respect to the muon. 



The effect of the first two factors is illustrated in figure pA\ . At 1 TeV the 
average difference between the u direction and the fi is about 0.7°. The difference 
decreases with increasing u energy. 




Figure 3.4: Angular difference between the initial neutrino direction and that of the 
muon track at the detector; the functional form shown reproduces the observed energy 
dependence well (solid curve). 



The detector resolution will be determined by the quality of the alignment of 
the detector components, the time resolution of the photomultipliers, the global 
timing of the readout system and the quality of the muon reconstruction. The 
reconstruction will be affected by light coming from secondary particles and by 
scattered light. Monte Carlo studies show that an angular resolution of 0.1° is 
possible. This means that above 100 TeV the total angular resolution is dominated 
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by detector effects whereas below 10 TeV the resolution is dominated by the 
angular distribution of the neutrino interactions. 



3.4.2 Energy response for interactions 

The energy response is determined by the energy fraction transferred to the muon 
in the neutrino interaction, the energy lost by the muon outside the detector and 
the energy resolution of the detector. The muon energy determination requires 
different techniques in different energy ranges. 

Below 100 GeV, the muons are close to minimum-ionizing, and the energy 
of contained events, with start and end points measured inside the detector, can 
be determined accurately from the range. The threshold for this method is about 
5-10 GeV for vertical tracks, depending on the vertical distance between groups 
of optical modules, and about 15 GeV for more isotropic events, depending on the 
horizontal distance between lines. 

Above 100 GeV, the range cannot be measured because of the limited size of 
the detector, but the visible range determines a minimum energy that can be used 
for the analysis of partially-contained events: starting events in which the vertex 
point is measured inside the detector, and stopping events in which the endpoint 
is measured. 

Above 1 TeV, stochastic processes (bremsstrahlung, pair production, (5-rays) 
are dominant, and the muon energy loss becomes proportional to the energy. The 
muon range above 1 TeV increases only logarithmically with the muon energy 
(figure |33| ). On the other hand, the detection efficiency increases with energy 
because of the additional energy loss. The correlation between measured muon 



energy and neutrino energy is shown in figure Monte Carlo studies have 
shown that the neutrino energy can be determined within a factor 3 above 1 TeV 
from the average energy loss. 

Above 1 PeV, the Earth becomes opaque to upward-going vertical neutrinos. 
Higher energies are accessible closer to the horizon, however. Very high-energy 
tau neutrinos can be observed because the produced in interactions decay 
before they are absorbed, producing u^. of lower energy which continue along the 
original z/^ flight path, but with decreasing interaction probability, resulting in an 
accumulation of events at the highest detectable energies. 
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Figure 3.5: Average range of the muon in standard rock as a function of the initial 
neutrino energy. 




Figure 3.6: Muon energy at the detector as a function of the parent neutrino energy. 



33 



3.5 Observable sky 



The ANTARES neutrino telescope, situated at a latitude of 43° North, can observe 
upward-going neutrinos from most of the sky (about S.Svr sr), due to the rotation 
of the Earth. Declinations below —47° are always visible, while those above 
+47° are never visible. Declinations between —47° and +47° are visible for part 
of the sidereal day (figure P77|). Most of the Galactic plane is visible, and the 




Figure 3.7: Visible sky in Galactic coordinates for a detector situated at 43° North 
latitude. The area within the dashed line is not observable; the area within the solid 
line is observable 24 hours a day. The sources shown are from the EGRET catalogue. 



Galactic centre is visible most of the sidereal day. Since the AMANDA telescope 
at the South pole is sensitive to positive declinations, the two detectors will have 
a reasonable area in common for cross-checks (about l.Svr sr). 

At energies greater than ~ 40 TeV, the interaction length becomes smaller 
than the Earth's diameter for z/^ traversing the dense core of the Earth. Above 
10 PeV, only nearly horizontal are visible (see figure If the field of view 
can be extended to 10° above the horizon at these energies where the background 
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is greatly diminished, a non-negligible fraction of the sky can be kept observable 
even at these energies. 




Figure 3.8: Transmittivity of the Earth as a function of incoming neutrino energy and 
zenith angle. 

Like other underground detectors, neutrino telescopes can observe the sky 
independently of the time of day, the phases of the moon or the weather. Existing 
underground experiments have usually reached 80% duty cycle after the initial 
debugging phase. ANTARES aims for even higher values for the off-shore 
facilities because the access is complicated and time-consuming. 
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Chapter 4 

R&D programme 



In order to ensure the success of the deployment of a large-scale detector in an 
uncontrollable environment such as the deep sea, it is necessary to perform an 
extensive programme of site evaluation and prototype testing. 

The selection of a suitable site for a neutrino telescope requires consideration 
of water transparency, optical background, fouling of optical surfaces, strength 
of the deep sea currents, meteorological conditions, depth, on-shore support, 
infrastructure and pier availability. This chapter reports on the results of site 
evaluation studies for the proposed site and on initial experience of prototype 
string construction and deployment. 

4.1 Site evaluation mooring lines 

A detailed programme of in situ measurements has been undertaken since 
October 1996, most of the data being taken at a site near Toulon (42°50' N, 
6° 10' E) at a depth of 2400 m (figure pTTj ). Optical background data were also 
taken 20 nautical miles off Porto, Corsica (42°22' N, 8° 15' E) at a depth of 
2700 m. 

Autonomous mooring lines to measure the site parameters have been 
developed and about 20 deployments and recoveries have been successfully 
performed. Three generic test lines exist, as illustrated in figure ^31 . Each test 
setup is incorporated in a mooring line anchored at the sea bed and vertically 
supported by a buoy. The measuring system, which can be up to 50 m long, 
is placed around 100 m above the sea bed. A further 30 m of cable and 
instrumentation sits between the measuring system and the top of the test line. 
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Figure 4.1: Map of the ANTARES site near Toulon. 

Tlie electronics, data acquisition and the detectors are powered by a set of lithium 
battery packs. A 500 cm'^ pack delivers 8 A-h at 26 V. These compact power 
sources have a lifetime which permits stand-alone tests for periods as long as one 
year. 

The test lines have been deployed using a standard horizontal deployment 
method. The buoy at the top of the line is immersed first, so that the line is 
held under tension during the entire operation. The boat continues to move slowly 
away from the buoy as the rest of the line is paid out from the boat. The anchor is 
the last element to be deployed after which the line falls to the sea bed. 

To date, the measurements carried out confirm that the properties of this site 
satisfy the constraints of the ANTARES physics programme. 

It is planned to explore a wider area in the Mediterranean sea, for long periods. 
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Figure 4.2: Test lines to study (a) light attenuation, (b) optical background, and (c) 
optical fouUng. 

in order to study the variation of parameters as a function of site, depth and 
season. For specific parameters such as water transparency, measurements at 
several wavelengths are also foreseen. The results obtained so far have been 
instrumental in steering the design of the telescope. 

4.2 Optical properties of the site 
4.2.1 Optical background 

The behaviour of the optical background on site places constraints on the trigger 
logic and the electronics as well as the mechanical layout of the optical modules. 
A setup has been devised for studying the time dependence of the background as 
well as its spatial extent and its correlation with deep sea current. The relevant 
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test line has been immersed ten times in total, for periods spanning from hours to 
months. 

The optical modules used for these tests are similar to those discussed in detail 
in section |5]3|. In this case 8-inch photomultiplier tubes were used. 

Up to three optical modules have been used on the same line, two of which, 
A and B, were 0.5 m to 1.5 m apart, while the third, C, was 10 m to 40 m 
away. A current-meter was installed below the optical modules. Data consisted 
of measurements of singles rates for all three optical modules and coincidence 
rates for modules A and B within a time window of 100 ns. In order to sample 
long-term variations, the system was enabled for a few hours three times a week. 
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Figure 4.3: Time dependence of the counting rate. 

An example of the observed counting rate is presented in figure ^j3| which 
exhibits two distinct components: a low level background around 40 kHz, and, 
superimposed on this, rapid (~ 1 s) excursions of up to several MHz. The low 
level background varies from 17 kHz to 47 kHz over a time scale of a few hours. 
This rate changes simultaneously on all optical modules even when they are 40 m 
apart. The peak activity is correlated with the current speed and is limited in 
spatial extent: peaks are seen simultaneously by optical modules when they are 
less than 1.5 m apart, but not when they are more than 20 m apart. The dependence 
of bioluminescence activity on current velocity is emphasized in figure where 
a correlation is observed between the two variables. This figure also illustrates 
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evidence for site to site and seasonal variations in this correlation. However, 
further data are required to understand these effects. 



BIOLUMINESCENCE vs. CURRENT SPEED 




current speed (cin/s) 

Figure 4.4: Bioluminescence burst activity as a function of current velocity. A parabolic 
fit is superimposed on the data. 



4.2.2 Optical fouling 

When exposed to sea water, the surfaces of optical modules are fouled by the 
combination of two processes: living organisms, mostly bacteria, grow on the 
outer surface, and sediments fall on the upward-looking surfaces. While the 
bacterial growth is expected to be almost transparent, sediments will adhere to it 
and make it gradually opaque, thus diminishing the sensitivity of the detector. This 
phenomenon is expected to be site-dependent as the bacterial growth decreases 
with depth and the sedimentation rate depends on local sources of sediments 
such as nearby rivers. A series of measurements has been performed in order 
to quantify these phenomena. 

Light of wavelength 470 nm from a blue LED source in a glass sphere was 
normally incident on a set of five PIN diodes placed at 50° to 90° from the vertical 
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axis of the sphere inside a second sphere 1 m away. Measurements of the light 
transmission and current velocity were made twice daily. 
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Time (days) 



Figure 4.5: Light transmission as a function of time and polar angle. 



Figure |4^ shows the light transmission monitored for 240 days with this 
configuration. For the horizontally-looking PIN diode (90°), the light source is 
affected by the same fouling as the detector, thus doubling the effect. In this case, 
a transmission loss of 1.2% per surface is observed after 8 months of exposure. 
All optical modules in the ANTARES detector will point downward, reducing the 
loss even further. 

In addition, a cradle holding glass slides mounted on a horizontal cylinder 
with various orientations around the axis of the cylinder (see figure was 
incorporated into the mooring line. This apparatus was developed by IFREMER 



Two series of measurements have been performed with the apparatus being 
immersed for long periods (3 months and 8 months respectively) at 2400 m depth 
on the same site as the other tests. After recovery, a biochemical analysis of 
the slides has been performed by several laboratories giving the results shown in 



figure 4-. 6. demonstrating the dependence of the biofouling on the orientation of 
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the slides. In those cases, where the slides point downward, saturation is observed. 
According to marine biologists, these numbers are several orders of magnitude 
below what is observed at shallow depths. 




Figure 4.6: (left) Position of glass slides on cylinder, (right) number of deposited bacteria 
vs. cylinder position 

The optical modules are housed in pressure-resistant glass spheres produced 
by the Benthos company. The mooring line which measured the fouling of 
Benthos spheres also included a sediment trap in order to determine the vertical 
flux and the composition of sediments at the ANTARES site. Particulate matter 
contributes to the scattering of light in sea water as well as the fouling of the 
Benthos spheres. The sediment trap collected samples from July to December 
1997 on a weekly basis. The analysis of the samples was performed in the 
CEFREM (Centre de Formation et de Recherche sur I'Environnement Marin) 



laboratory by the team of Professor A. Monaco. As shown in figure |4J| , the total 
mass flux of sediments substantially increases from October onwards, probably 
as a result of heavy rainfall draining sediment from the shore. Indeed, the 
composition of the sediments shows a large contribution of material originating 
from continental river bed. This study will be complemented by the analysis of 
the sediment cores and of the water samples collected during the ANTARES site 
survey campaign of December 1998. 
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Figure 4.7: Total mass flux of sediments at the ANTARES site as a function of time from 
July to December 1997. 

4.2.3 Transmission properties of the water 

The water transparency affects the muon detection efficiency, while the amount of 
scattered light determines the limit on the angular resolution of the detector. These 
two parameters therefore influence the detector design and are necessary data for 
the Monte Carlo simulations from which the detector response is calculated. They 
must be measured in situ, as water samples may be degraded when brought to 
the surface. Two different experimental setups have been constructed to measure 
these parameters. 

In December 1997, measurements were performed with a 33 m long rigid 
structure holding a coUimated and continuous LED source located at a variable 
distance from an optical module. For each selected distance D between the source 
and the detector, the LED luminosity $led was adjusted so as to yield a constant 
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current /pmt on the photomultiplier tube. The set-up was calibrated with a similar 
experiment done in air. The emitted and detected intensities in water being related 
by 



JpMT oc <I>led/-D X exp(-Z:'/Aatt.cfr) 



(4.1) 



this test makes it possible to estimate the effective attenuation length from the 
dependence of the required LED intensity with the distance (cf. figure Ol). The 
agreement of the data with a decrease following the formula given above was 
excellent and yielded an effective attenuation length of 



Aatt.off = 41 ± 1 (stat.) ± 1 (syst.) m (December 1997) 



(4.2) 



This attenuation length results from a combination of absorption and scattering. 
The experimental set-up was unable to separate these, and the long rail made 
deployment difficult. The experiment was therefore redesigned so that it used a 
pulsed source, to facilitate scattering measurements, and a flexible structure to 
improve ease of deployment. 




• Measurement in water 
Calibration in air 
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Figure 4.8: Determination of attenuation length. The log of D^/flux is plotted against D, 
the distance between the LED and the optical module. 

Li July 1998 and March 1999, measurements were performed with a set-up 
consisting of a pulsed isotropic LED source located at a distance of either 24 m or 
44 m from a 1" fast photomultiplier tube. An 8-bit TDC measured the distribution 
of the arrival times of the photons. The overall time resolution was a = 4.5 ns. 
Because of this, photons in the tail of the distribution have scattered with an angle 
at least ~ 35° for the 24 m spectrum, and ~ 25° for 45 m. Therefore, the scattering 
properties of the water are being measured for large scattering angles. 
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Figure 4.9: (left) Distribution of arrival times of photons for two distances between the 
detector and the source. The 24 m distribution is normalized to 1 , and the 44 m distribution 
is normalized relative to the first one. (right) Arrival time distributions for 24 m, 44 m, 
and a calibration in air, with Monte Carlo curves superimposed. 

The time distributions recorded exhibit a peak stemming from direct photons, 
and a tail extending to larger delays due to scattered photons. For the 24 m (44 m) 
spectrum, 95% (90%) of the photons are collected within 10 ns. Scattering is thus 
a small effect at the ANTARES site, as illustrated in figure P~9| . 

An effective attenuation length could be determined from the ratio of the 
integrated spectra measured at the two distances, yielding: 

60.0 ± 0.4 (stat.) m (July 1998) 
Aatt,efr- 52.2 ±0.7 (stat.) m (March 1999) ^ 

A systematic uncertainty of a few metres might affect these estimates due to the 
fact that the LED luminosity is not monitored and yet assumed to be the same for 
the time distributions collected at the two distances. These measurements indicate 
more significant attenuation in March than in July. The difference between these 
and the December 1997 measurement is also significant, although it is partly 
accounted for by the use of a coUimated source in 1997. 

A more detailed analysis can be done by fitting the data with a Monte Carlo 
distribution obtained by photon tracking. The data are well described using an 
absorption length in the range 55-65 m, a scattering length at large angles greater 
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than 200 m and a roughly isotropic scattering angle distribution. This is consistent 
with the effective attenuation length deduced from the ratio of integrated spectra. 

4.3 Sea conditions 

Suitable sea conditions for periods of up to a few consecutive days are required 
to perform deployment and recovery operations. These conditions depend both 
on the nature of the operations and on the characteristics of the ship. For the 
single string deployment and recovery operations of June-September 1998 with 
the Castor, a wave height less than 1.5 m and wind speed less than 25 knots (5 on 
the Beaufort scale) were specified. 

A study has been made incorporating data from a number of sources, namely: 

• data on wave height collected by an instrumented buoy moored 4 nautical 
miles south of PorqueroUes Island from May 1992 to September 1995, 
which should be representative of the conditions in the ANTARES site; 

• data on the wind speed and direction as recorded by the PorqueroUes Island 
Signal Station analysed for the period from January 1992 to March 1996; 

• additional information on 
sea conditions provided by satellite measurements in the area from 1992 
onwards. 

Preliminary analysis of these data has been performed by the Meteomer 
company. Periods of three consecutive days with favourable sea conditions occur 
less than five times per month from October to April, and more than five times per 
month from May to September. 



4.4 Site survey 

The strength and direction of the underwater currents need to be taken into 



account in the mechanical design of the detector. Figure [4-.10| summarises all the 
measurements of the deep sea current gathered during the test immersions. The 
maximum current observed, namely 18 cm/s, is accommodated in the mechanical 
design of the strings. 

A visual and bathymetric survey of the sea floor was performed in December 
1998, using the Nautile submarine. In the area selected as a potential ANTARES 
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Figure 4.10: Summary of the survey of deep sea currents. 

site, the sea floor is flat and displays no topographical anomalies such as steps 
or rocks, as can be seen in figure [4.1 1| . During the same series of dives, core 
samples from the sea floor were collected. They consist of solid mud which is a 
satisfactory substrate to support the detector. 

4.5 Prototype string 

To gain realistic experience of the issues involved in the deployment, operation 
and recovery of a full-scale detector, a full-size prototype of an ANTARES 
detector string has been constructed. It is equipped with a positioning system, 
slow control network, power distribution, and eight optical modules with their 
associated readout electronics. Once deployed, the string is foreseen to be linked 
to a shore station by an electro-optical cable supplying the power to the string as 
well as enabling the control and readout connection. The string is adequate for a 
detailed study of the optical background and a measurement of the down-going 
muon flux. 

4.5.1 Mechanical structure 

The prototype string is 350 m high (figure [4-.12[ ), anchored on the sea floor and 
vertically supported by a buoy. It is composed of two vertical cables spaced 
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Figure 4.11: Site survey map 



2.3 m apart, supporting both ends of sixteen Optical Module Frames (OMFs). 
The OMFs are placed every 15m from a height of 95 m up to 320 m from the sea 
floor and are constructed from fibre glass to avoid corrosion. Horizontal spacers 
keep the vertical support cables below the lowest OMF; a spacer is also situated 
at the middle of every 15 m segment separating the OMFs. 

Each OMF supports a pair of Optical Modules, separated by 1.6 m from axis 
to axis, while the central part houses a container for electronics made of corrosion- 
resistant titanium alloy. 



4.5.2 Slow control system 

A slow control network is linked to the shore station through an electro-optical 
cable and employs an architecture whereby each electronics container has a point- 
to-point connection with the slow control data acquisition system. It permits 
the control of the power distribution electronics, optical module motherboard. 



49 



650 daN buoyancy 



10 01 Optical module frames 
# Optical modules 

Instrumented spheres 

™ satellites 

□ LCM 

MEC 

1 hydrophones 

□ celerimeter 

H acoustic container 
CTD 

I^H currentmeter 



Electro-optical cable 



[05Q] 



fo^ 



4^ ^^J ^ 





340 m 



— 306.6 m 



15 m 



102.2 m 
87.6 m 



m 



2000 daN anchor 



Figure 4.12: Schematic view of the prototype string. 
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analogue readout electronics and the acoustic positioning system. It also handles 
the readout and transmission of sensor data, such as the satellites, the acoustic 
system and the electronic card temperatures or status. 

This DAQ system is located in the main electronics container at the bottom 
of the string. The main fibre-optic data link connects this system to the shore 
station, which receives the slow-control data and provides a user interface for the 
slow-control system. 



4.5.3 Positioning 

A string does not provide rigid support for the optical modules. Two independent 
systems have been incorporated in the prototype string to provide a precise 
knowledge of the relative position of each OM at any time. The first system is 
based on a set of tiltmeters and compasses which measure local tilt angles and 
orientations on the string. The reconstruction of the line shape, as distorted by 
the water current flow, is obtained from a fit of measurements taken at different 
points along the line. A successful test of this system was performed during a 
deployment of the prototype string which was equipped with several sensors with 
a precision of 0.05° in tilt and of 0.3° in direction. Figure presents a view of 
the reconstructed line, the segments indicating the sensor positions and directions 
with respect to North. A maximum error of ~ 1 m on the reconstructed shape is 
estimated. 

The second system, based on acoustic triangulation, is more precise but 
requires more complex and expensive electronics. In this system, rangemeters 
placed on the string send an acoustic signal to a minimum of three transponders 
fixed to the sea bed. Each transponder replies with its characteristic frequency. 
A global fit of the measured acoustic paths gives the precise three-dimensional 
position of the rangemeters, provided that the positions of the transponders and 
the sound velocity in water are known. The prototype string is equipped with 
four rangemeters (a hydrophone with its electronics container) communicating 
with four external autonomous transponders placed on a fixed structure on the sea 
floor about 200 m away from the string. The measurement of the communication 
time between one fixed rangemeter and one fixed transponder demonstrates a 



reproducibility of ~ 1 cm in the acoustic path length (figure |4.14j ). 

In order to exploit such a system fully, a precise knowledge of the sound 
velocity in water along the acoustic path is required. This depends strongly on 
water temperature and also on salinity and depth. The prototype string is thus 
equipped with sound velocimeters, which measure the local sound velocity with a 
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Figure 4.13: Top view of the reconstructed shape of the prototype string during 
deployment, sensors are labelled sat 1 to sat 9 

precision of 5 cm s~^, and with Conductivity Temperature Depth devices (CTDs) 
to observe the variations of temperature and salinity. 

The systems are complementary: a few points of the line can be measured 
acoustically and other points are obtained by line shape fitting. The tests already 
performed confirm that the desired precision on relative OM positioning can be 
achieved. 

4.5.4 String deployment 

The prototype string was deployed, operated and retrieved several times in 
Summer 1998, first at 400 m depth and then at 2300 m. This string was powered 
by batteries housed in a container at the bottom of the string. The deployment 
procedure used was different from that used for the mooring lines discussed 
earlier. In this case, a step-by-step procedure involving two winches on the 
boat deck was implemented whereby the anchor was immersed first, the string 
being held securely at 2 points, then each OMF was paid out storey-by-storey. 
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Figure 4. 14: Distribution of the acoustic path measurement between a fixed rangemeter 
and a fixed transponder. 

Recovery was performed in a similar way. Deployment and recovery at 2300 m 
took 18 hours in total. This method allowed the equipment to be deployed in a safe 
and controlled manner. The exercise also permitted the handling procedures for 
the electro-optical cable termination attached to the string anchor to be verified. 
The tiltmeter and compass data were recorded so that the detector string behaviour 
could be studied during the deployment phase. The fully equipped version of the 
prototype string will be ready for immersion during Summer 1999. 

Successful deployment of the 0.1 km^ array will require strings to be located 
in well-defined positions with 80 m spacing or less. This was investigated in a 
multi-string deployment test using the dynamical positioning ship Provence and 
400 m test strings instrumented with acoustic beacons. 

String positioning precision was investigated by aiming to deploy a second 
test string 50 m away from the one which had been deployed earlier. Relative 
positions were measured with the aid of three external acoustic beacons. 

The second string was deployed using a winch, with the ship located at the 
target position. Deployment was halted with the string 200 m above the sea floor. 
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the relative position measured, and the ship moved to adjust as necessary. Using 
this technique it proved possible to locate the second string to within 10 m of the 
nominal position. 



4.5.5 Deep sea connection and line recovery 

In the full-scale detector, each string will be connected to a common point, known 
as the junction box. A connection scheme from the junction box to a prototype 
string anchor was tested during a ten-day sea operation at the ANTARES site, 
using the Nautile submarine and its support vessel the Nadir from IFREMER. 

A reel of cable equipped with deep sea connectors on both ends was immersed 
first. Then the Nautile uncoiled the cable by pulling it and plugged the connector 
of the cable to its counterpart on the anchor. The connection procedure was 
successfully performed twice, using a free-flight technique which requires only 
one of the submarine's arms. This offers the possibility of using IFREMER's 
second manned submarine, the Cyana, which is equipped with only one arm. The 
advantage of this is that this submarine is more freely available than the Nautile. 



Figure |4-.15| shows a diagram of the Nautile approaching the anchor, holding 
the connector and the cable, before plugging it at the end of the arm which can 
be seen in an upright position at the bottom of the string. Once the connection 
is made, the Nautile pushes the arm down and plugs it into a fork linked to the 
anchoring weight. To retrieve the string, acoustic releases are activated from the 
surface which disconnect the string from the anchoring weight. The buoyancy 
of the string pulls on the connector which is held back by the fork and is thus 
unplugged. This system avoids the need for a submarine for string retrieval. 

The speed of ascent of the string during its trip to the surface is around 1 m s~^. 
The distribution of weights and buoyancies along the string have to be carefully 
studied in order to control the relative speed of each storey so that the string does 
not become entangled when it surfaces. 
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Chapter 5 
Detector design 



The detector design described in this chapter is constrained by the following 
factors: 

• the environmental conditions at the site, which influence the spacing of 
optical modules, the mechanical structure of detector strings, and the 
expected background rates; 

• the practical experience gained from the design and deployment of the 
protoype string; 

• the need to optimise the physics performance of the detector; 

• the requirement for a high level of reliability. 

The first two issues were discussed in the preceding chapter, and the 
physics performance and its optimisation are the subject of chapter ^. A key 
issue throughout the detector design is the question of reliability. Repair and 
maintenance of offshore detector elements will involve high costs and probably 
long down times, and must be minimised. For this reason, the design of each of the 
detector components discussed in the subsequent sections will follow a strategy 
intended to enhance reliability: reduction of the number of active components, 
limitation of power consumption, and avoidance of single-point failure modes. 

This chapter outlines the specifications satisfying these constraints. The 
design may evolve further as technical details are refined and additional physics 
simulations are performed. 
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Figure 5.1: Schematic of part of the detector array; the magnified view shows two storeys 
and a hydrophone. 



5.1 Overview 



The detector consists of an array of approximately 1000 photomultiplier tubes in 
13 vertical strings, spread over an area of about 0.1 km^ and with an active height 
of about 0.3 km. Figure ^]T] shows a schematic view of part of the detector array 
indicating the principal components of the detector. 

The basic unit of the detector is the optical module, consisting of a 
photomultiplier tube, various sensors, and the associated electronics, housed in 
a pressure-resistant glass sphere. The electronics includes a custom-built digital 
electronic circuit which captures and stores waveforms, pulse heights and timing 
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information, as well as the HV power supply for the photomultiplier tubes and the 
network nodes for data transmission and slow control. 

The optical modules are grouped together in 'storeys' of three modules and 
interconnected via an electro-mechanical cable. In the present design the detector 
has 13 strings, each of which has a total height of about 400 m. Of the 13 strings, 
four consist of 41 storeys spaced vertically by 8 m and nine have 21 storeys spaced 
vertically by 16 m. The strings are arranged in a randomized spiral (see figure |677|) 
with a minimum horizontal spacing fixed at some value between 60 m and 80 m. 
The smaller spacing gives a higher detection efficiency for muon energies up to 
50 TeV but could be more difficult to deploy. 

The optical modules in a storey are arranged with the axis of the 
photomultiplier tubes 45° below the horizontal. The present plan is to use 10- 
inch Hamamatsu photomultiplier tubes, but other, possibly larger, tubes are under 
development by various manufacturers. The angular acceptance of the optical 
modules is broad, falling to half maximum at ±70° from the axis (see below). 
This means that the proposed arrangement of OMs detects light in the lower 
hemisphere with high efficiency, and has some acceptance for muon directions 
above the horizontal. In the lower hemisphere there is an overlap in angular 
acceptance between modules, permitting an event trigger based on coincidences 
from this overlap. 

The relative positions of all optical modules in the detector are given in real 
time by a positioning system identical to that described for the prototype string in 
section [4.5.3 . 

Each string is instrumented with several electronics containers. At every 
storey, there is a local control module (LCM), and at the base of each string there is 
a string control module (SCM). Special containers house acoustics and calibration 
equipment. Each of these containers constitutes a node of the data transmission 
network, receiving and transmitting data and slow-control commands. The 
functions which they support include reading sensors, adjusting slow-control 
parameters, the trigger, and the distribution of power, master clock and reset 
signals to the front-end electronics. 

The individual SCMs are linked to a common junction box by electro-optical 
cables which are connected using a manned submarine. A standard deep sea 
telecommunication cable links the junction box with a shore station where the 
data are filtered and recorded. 

The trigger logic in the sea is planned to be as simple and flexible as possible. 
The first-level trigger requires a coincidence between any two OMs in a single 
storey. The second-level trigger is based on combinations of first-level triggers. 
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Following a second-level trigger the full detector will be read out. A more refined 
third-level trigger, imposing tighter time coincidences over larger numbers of 
optical modules, will be made in a farm of processors on shore. The readout 
rate is expected to be several kHz, and the corresponding data recording rate less 
than 100 events per second. 

The following sections of this chapter describe the various components of the 
detector in more detail. 

5.2 Detector string 

The design chosen for the detector string is similar to that used for the prototype 
string described earlier in chapter ^, a string maintained vertically by its own 
buoyancy and anchored on the bottom of the sea. Between the buoy and the 
anchor, the active detector part of the string comprises a series of elementary 
detector segments. These segments are standardised and so can be mass-produced, 
and if necessary interchanged. 

In all cases, the design of the string components meets basic specifications 
such as corrosion resistance, required because of the prolonged immersion in salt 
water, resistance to high pressure and water-tightness. Similarly, in certain cases, 
the colour and surface properties of detector elements need to be selected so as to 
minimize light reflection. All detector components must remain functional for a 
minimum lifetime of 10 years. 

The principal elements of each string are illustrated in figure ^]2| and discussed 
further here. 

The bottom string socket (BSS) anchors the string to the sea bed, facilitates 
the electrical connection of the string to the network and permits the release 
and subsequent retrieval of the string. To allow precise and simple string 
installation its construction has been optimised for handling on the deployment 
ship, resistance to shock, stability during descent. The BSS is instrumented for 
acoustic positioning. Connection of the string to the network is performed by a 
submarine and is described in section [4-.5.5| . 

The electro-mechanical cable (EMC) provides mechanical support for the 
string as well as enabling the electrical interconnection of the detector string 
elements. It must be capable of supporting tensile, torsion and bending stresses in 
order to maintain the string's stability. Its construction must be suitably flexible 
to allow the integration of the various configurations, as well as supporting the 
service, storage, handling and immersion and recuperation phases. Electrical 
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Figure 5.2: Sketch of the detector string, abbreviations are defined in the text. 
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cables and optical fibres run through the EMC. They enable power distribution and 
the transmission of signals between two consecutive electronic containers (LCM 
or SCM). 

Each elementary segment (ES) consists of three optical modules, one local 
control module and one optical module frame. The OMF supports the various 
elements placed within it, namely the optical modules, the LCM container and 
possibly a hydrophone for acoustic positioning. The OMF design must be 
sufficiently flexible to accept certain changes in the number and positioning of 
the optical modules. Major traction forces that work on the substructure during 
deployment and recovery must in no way be transmitted to those string elements 
supported by the frame. 

The top of the string consists of a buoy. Its dimensions and geometry are 
optimised to minimise hydrodynamic effects such as dragging and vibration whilst 
maintaining a suitable tension in the string. In addition, the buoyancy should be 
sufficient to ensure a controlled resurfacing of the string during retrieval. 



5.3 Optical module 



A schematic view of the ANTARES optical module is shown in figure ^ 
Mechanical and optical aspects of the OMs are reviewed in this section, as well as 
the performance of the various photomultiplier tubes which have been evaluated. 



5.3.1 Benthos sphere 

The photomultiplier tube and its associated electronics are housed in a 43 cm 
diameter, 15 mm thick. Benthos sphere, which can withstand pressures of up to 
700 bars. The sphere is made of two halves, one of which is painted black on its 
inner surface so as to give the OM some minimal directionality with respect to 
Cherenkov light detection without degrading its acceptance. The two halves of 
the sphere have machined edges which form a seal when subjected to an external 
over-pressure. Attenuation of light at A = 450 nm due to the sphere was measured 



to be less than 2% (see figure 5.4). 



5.3.2 Optical contact and magnetic shielding 

Silicone gel ensures both optical coupling and mechanical support of the PMT. 
The refractive index of the gel (Ugei = 1.40) does not exactly match that of 
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Figure 5.3: Cross section of the optical module. 



the sphere itself (n glass = 1.48) but is higher than the refractive index of water 
(neater = 1-35) and so the amount of light reflected out of the OM is minimized. 
The silicone gel covers the entire photocathode area; its attenuation length is given 
as a function of wavelength in figure ^]4| [|7]] . 

The Earth's magnetic field significantly degrades the collection efficiency 
of phototubes by bending the trajectories of electrons, mainly between the 
photocathode and the first dynode. A cage made of 1 . 1 mm thick high-permittivity 
alloy wire is used to shield the PMT and to minimise the dependence of the OM 
response with respect to its angle with the magnetic North. The mesh size of the 
cage (6.8 cm) was optimized to reduce the non- uniformity of the PMT angular 
response to less than 5% while minimizing the fraction of light lost due to the 
shadow on the photocathode 



5.3.3 Photomultiplier tube characterisation 

During the last three years several test benches have been set up to measure 
and compare the principal features of different PMTs and to provide detailed 
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characteristics of the OM response. These include several 'dark boxes', in which 
the PMTs are exposed to uniform illumination coming from red, green or blue 
LEDs or very fast solid-state lasers. This allows the systematic and precise 
measurement of a large number of PMT characteristics [|60|]. Several PMTs of 
different diameters have been studied, and their performances are summarised 
here in terms of a small number of critical parameters. This discussion only 
considers a fraction of the total number of parameters which have actually been 
measured (see [g^, ^, ^). 

The effective photocathode area (^^/) is defined as the detection area of the 
photocathode weighted by the collection efficiency. It was measured by scanning 
the entire photocathode surface with a coUimated blue LED. 




Figure 5.4: Attenuation length of (a) the Benthos sphere, (b) silicone gel and (c) sea 
water; (d) quantum efficiency for the Hamamatsu photomultiplier tubes. 



Electromagnetic interference in the optical module induces noise at the PMT 
anode. It is expected that this will not exceed 5 mV (rms). This noise governs 
the gain at which the PMT is operated. A factor of 10 between the average 
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pulse height for a single photo electron (SPE) and the noise is sufficient to ensure 
efficient discrimination of the signal. This corresponds to an effective working 
gain of the order of 5x10''. In view of PMT ageing and possible variations in 
specifications a maximum gain of at least 10^ is required. 

The Peak to Valley ratio is computed from the observed charge spectrum of 
single photoelectrons with the high voltage adjusted to give 50 mV amplitude for 
SPE. The Peak to Valley ratio is required to be greater than 2. 

Due to imperfections in the electron optics and the finite size of the 
photocathode, the SPE transit time between the photocathode and the first dynode 
has a measurable width, usually referred to as the transit time spread (TTS). This 
defines the timing resolution of the PMT, which is required to be comparable to 
that from the overall positional accuracy and the timing precision in the readout 
electronics, i.e. 1.3 ns rms or 3 ns FWHM. The measurement of the TTS is 
performed over the whole photocathode area with the PMT operating at a gain 
of 5x10^. 

The maximum quantum 



Four different PMTs are compared in table 



efficiency, the linearity and the dark count rate have also been determined. 



Manufacturer 


Hamamatsu 


ETL 


Phototube diameter 


8" 


10" 


8" 


11" 


^f/'/ (cm^) 


280 


440 


240 


620 


Maximum gain 


109 


109 


108 


5x10^ 


Peak to valley ratio 


2-3 


3-3.5 


< 2 


2 


TTS (FWHM) (ns) 


2.5 


3.5 


2.5 


3.0 



Table 5.1: Comparison of relevant properties of large-photocathode phototubes tested to 
date. 



5.3.4 Optical module response 

To reproduce the experimental conditions, it is necessary to investigate the 
response of an optical module to the Cherenkov light produced by muons 
travelling through water. This has been studied by placing an OM in a 1 .5 m high, 
1.4 m diameter, light-tight cylindrical steel tank filled with constantly recycled 
fresh water [^. Nearly- vertical atmospheric muons are tracked in a four plane 
hodoscope, two planes above the tank, two below. Between the tank and the two 
lower planes a lead shield ensures Ef^ > 0.6 GeV Muons traverse the tank and 
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produce Cherenkov light which illuminates the optical module. The light yield 
for muons of this energy is more than 90% of the value for a single relativistic 
track. The optical module is able to rotate, so that its response can be measured at 
different angles, Oqm, between the incident muon and the PMT axis. 

Figure |5]^ presents results from this setup using an optical module equipped 
with a 10" tube. The angular region where the OM sees more than half the 
maximum amplitude is of the order of ±70°. 




20 40 60 80 100 120 140 160 180 



Sqm (degrees) 



Figure 5.5: Number of photo-electrons detected by the optical module as a function of 
its zenith angle Oqm, for nearly vertical muons normalised at 1 m distance. 



5.4 Offshore electronics 

The distance between the detector and land precludes a point-to-point connection 
between each optical module and the shore station. Instead, an electro-optical 
cable from the shore station supplies electrical power to the detector array 
and permits data to flow in both directions. The electro-optical cable ends at 
the junction box to which the strings are connected. A star-topology network 
architecture is used, running from the string control module to the optical modules 
via the local control modules. A digital scheme has been developed for the 
necessary data multiplexing. 
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This network is used to distribute power, collect data, broadcast slow control 
commands and master clock signals, and form the trigger. This section describes 
the implementation of the above scheme. 

5.4.1 Digital front end 

The OM electronics must fit in the limited space available in the OM, consume 
little power, be reliable and long-lived (average lifetime > 10 yr), be inexpensive. 

An Application Specific Integrated Circuit (ASIC) meets these requirements 
and can be tailored to our needs. The ASIC developed for the digital front end 
is called the Analogue Ring Sampler (ARS), similar in design to the Analogue 
Transient Waveform Recorder []62|]. The chip samples the photomultiplier tube 
signal continuously at a tunable frequency between 300 and 1000 MHz and holds 
the analogue information on 128 switched capacitors when a low-level threshold 
is crossed. The information is then digitized by an external 8-bit ADC. Figure [5^ 
shows the resulting histogram for an effective sampling frequency of 1000 MHz. 
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Figure 5.6: Charge measured in the 128 switched capacitors of the ARS memory after 
the capture of an analogue signal from a photomultiplier tube. The reflection of the signal 
comes from an impedance mismatch on the test board. 

A 20 MHz reference clock is sampled on one channel, giving a relative timing 
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of the signals to better than 1 ns. A time stamp is obtained for each event by 
counting the reference clock cycles. A reset command sent through the clock 
stream is used to restart all the counters of the array synchronously. 

Since 99% of the pulses are single photo-electrons, another ASIC, dedicated 
to the treatment of single photo-electron pulses, has been developed in order to 
reduce the dead time and the data flow. The first part of the ASIC performs pulse 
shape discrimination (PSD). It identifies three types of pulses requiring analysis 
of the full wave form: 

• large pulses, with pulse heights that cross a threshold corresponding to 
several photo-electrons; 

• wide pulses, with time over a low threshold longer than about 15 ns; 

• two pulses separated by less than 50 ns. 

If none of these conditions is met, only the pulse charge and time of arrival are 
measured, this information, along with the OM address, is transmitted to shore 
in 64 data bits. The large pulse threshold, the time over threshold, and the time 
window for multiple pulses are adjustable. The remaining 1% of the pulses satisfy 
one of the three conditions above and so the pulse shape is transferred for offline 
analysis, approximately 2000 data bits are required to encode this type of event 
(waveform event). 

A new version of the ARS is currently under development. It integrates all 
these functions on the same chip, together with the ADCs, the DACs and the 
slow-control interface (figure p77| ). A pipeline memory is implemented to store 
the single photo-electron information long enough to match the level 2 trigger 
propagation and formation time, which is around 10 ;us for a 0.1 km^ detector. 
It will be possible to use up to four memories per optical module, in a token 
ring. This will permit the chip to be used with photomultiplier tubes bigger than 
10 inches in diameter, for which the counting rate may exceed 60 kHz, and for 
a km-scale detector, where the trigger formation and propagation time may reach 
30 fis. 

5.4.2 Trigger logic and rates 

A level trigger occurs when the output of a photomultiplier tube crosses a 
threshold corresponding to 30 % of a single photo-electron amplitude. In order to 
deal with the high counting rate in the sea (see section [4.2. 1[ ), a level 1 trigger is 



68 



Time ref. . 
clock 

Dyn. 2 . 
Dyn. 1. 
Anode . 



Waveform mode control (ARSO) 



Wavefomi mode 
- 4*128 memory cells - 
(ARSO) 



SPE mode 
Charge Integrator 



PSD 
manager & 
LO trigger 



SPE mode 
TVC 



SPE mode Pipeline 

WW 

— Analog memory cells — 



Pipeline control 



ADC manager 



ADC 



ADC 



DACs 
regisiers. 

configuration 
interface. 

Slow control 



•- Serial 
link 



► Serial 
link 



triggers 
LI, L2 
manager 



. LI trigger 
. L2 trigger 

Time ref. 
clock & 
Synchro 



Figure 5.7: ARSl block diagram. 



built out of a tight time coincidence between two level triggers from the same 
storey. A level 2 trigger can be formed by requiring multiple level 1 triggers in 
a coincidence gate whose width is of the order of that needed for a track to pass 
through the entire detector. The level 2 trigger condition could be at least two 
level 1 triggers on the same string (referred to as the 'string trigger') or at least 
three level 1 triggers anywhere in the detector (referred to as the 'array trigger'). 
When either the string trigger or the array trigger conditions is satisfied, a readout 
request is sent to the entire array. The readout request is received in each OM, 
which starts the digitization of all the information within the maximum allowed 
time of flight. Internal delays specific to each OM compensate for the trigger 
formation and readout request propagation times. When a level 1 trigger occurs 
in a storey, the two OMs involved are read out, even if they do not participate in 
the subsequent level 2 trigger. 

The level 1 trigger logic will be installed in each LCM. The level 2 trigger logic 
must be linked to all the LCMs which may participate in the trigger: therefore the 
array trigger will be installed in the junction box, while the string trigger will 
be installed at the bottom of each string, and then sent to the junction box, from 
which the readout request will originate. 

Trigger rates due to random coincidences from background counting rates 
have been estimated assuming a level trigger rate of 60 kHz. The level 1 trigger 
requires a coincidence between two of the three optical modules on any storey 
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with a pulse width of 20 ns, yielding a rate about 500 Hz per storey. 

Level 2 trigger rates are calculated assuming four dense strings of 41 storeys 
with 8 m vertical spacing, and nine sparse strings of 21 storeys with 16m spacing. 
The gate width for string triggers is 1.5 /iS, the time for a particle to travel 
the length of the string. For array triggers, which involve multiple strings, a 
coincidence width of 2 /iS is used. A simple string trigger requiring two level 1 
triggers on the same string leads to a level 2 trigger rate of about 3 kHz. A similar 
array trigger requiring three level 1 triggers within the 353 storeys of the full 
detector adds 6 kHz, yielding a total level 2 trigger rate of about 10 kHz. 

More restrictive level 2 triggers can be obtained by combining signals from 
separate sectors on a single string to yield 'proximity' triggers, favouring nearby 
hits on the same string. The time for light to traverse two sectors is 400 ns. 
Requiring two level 1 triggers within two adjacent sectors leads to a total rate 
of 300 Hz, down a factor 10 from the simple string trigger described above. 

Even more restrictive level 2 triggers can be defined if necessary. The easiest 
way to reduce the trigger rate is to increase the number of level 1 triggers required 
within the level 2 coincidence gate. For the simple string triggers, requiring three 
level 1 triggers would reduce their contribution to the level 2 trigger to about 
30 Hz, a reduction of a factor 100. For the array triggers, requiring four level 1 
triggers in coincidence would reduce the level 2 rate to about 600 Hz, a reduction 
of a factor 10. 

During a bioluminescence burst, level 1 triggers from affected storeys provide 
no discrimination against noise, because all the OMs in the storey see the burst. 
These storeys are removed from the trigger logic in real time, but the OMs are still 
read out, as the timing resolution of the detector is sufficiently precise to allow real 
hits to be recovered from the noise when tracks are reconstructed. 

5.4.3 Data flow rate 

The volume of data transmitted to shore depends on the trigger rate, the OM 
background and the proportion of waveform events. The following calculation has 
been performed for a total trigger rate of 10 kHz with a 2 /iS time window, so that 
2% of the overall background activity is read out and sent to shore. The noise rate 
from each OM is assumed to be 60 kHz for 95% of the time (quiescent phase) and 
300 kHz for the remaining 5% (bioluminescence activity). It is also assumed that 
the level 1 trigger is disabled for modules undergoing bioluminescence activity. 

The three classes of waveform events, namely large pulses, wide pulses and 
events where two pulses occur within 50 ns are discussed in section |5.4.1| . Large 
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or wide pulses originate from decays in or near the OM and amount to a few 
hundred Hz, while closely spaced pulses can result either from pre- or after- 
pulsing or from accidental coincidences. At 60 kHz, the three contributions are 
of approximately equal importance and give about 1 kHz of waveform events, 
while at 300 kHz, accidental coincidences dominate, yielding about 10 kHz of 
waveforms. 

Given these figures, the data flow rate from the level 2 triggers is 300 Mb/s. 
Reading out the OMs involved in level 1 triggers adds a further 100 Mb/s to give 
a total data flow rate of about 400 Mb/s. 

5.4.4 Data handling and transmission to the shore 

A number of electronic functions must be developed in order to handle the digital 
data coming from the optical modules and transmit it to shore. Functions are 
required at each storey on the string (LCM), at the base of the string (SCM), 
and in the junction box connected to the shore by the electro-optical cable. The 
specific functions are: 

• Organizing data coming from the optical modules and transferring them to 
the SCM. 

• Processing the level 1 trigger at the LCM and the level 2 trigger at the SCM 
and at the junction box. 

• Slow control and monitoring hardware. 

• Clock distribution from the shore to the LCM containers and on to the 
optical modules, in order to time stamp all digitized physical events. 

• Handling of the acoustic positioning electronics for each string, distributed 
over selected storeys throughout the length of the string, performed at the 
LCM level. 

• Acquisition of mechanical positioning information from tilt-meters and 
compasses incorporated in the LCMs. 

• Transmission of data and slow-control commands from the shore to the 
SCMs via the junction box. 

Experience with the prototype string has led to the adoption of the following 
techniques: 
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• Fibre-optics along and between strings, providing reliable, compact and 
high bandwidth links 

• Penetrators as opposed to connectors at the level of each LCM: this 
permits the use of optical fibres, allows connections to be made inside the 
container, and increases the number of communications channels which can 
be integrated into a single electro-mechanical cable. 

5.4.5 Electro-optical cable and junction box 

Existing electro-optical cables designed for submarine telecommunications 
applications fulfill the power and data handling requirements summarised above. 
They incorporate a copper coaxial conductor for power transmission and 16 to 24 
optical fibres for digital signals. 

The junction box contains power converters in an oil bath to provide a standard 
400 V DC voltage, the trigger electronics and the electro-optical interface for the 
data, slow control and clock transmission. 

The cable and junction box are a potential single point failure of the detector. 
It is therefore foreseen to have two cables and two junction boxes, each capable 
of handling the full requirements of the detector. 

Electro-optical cables link the junction box to each SCM. since optical fibres 
are necessary for the data transmission and the clock distribution. In order to use 
standard undersea electrical connectors, it is foreseen to convert all optical signals 
to electrical ones near the connectors. 

5.5 Slow control and commands 

The system of slow control is intended for the monitoring of variables which 
change relatively slowly, and also serves to control various aspects of detector 
operation. PMT voltage, temperature and power-supply voltages are read from 
the optical modules. Dedicated instruments provide information on string attitude 
and orientation, water current velocity, acoustic positioning information, and other 
control data. Parameters to be adjusted during detector operation include the 
PMT voltage, thresholds involved in pulse detection and triggering, and various 
calibration systems. The slow control system thus gives a user at the shore station 
all the information needed to monitor and control the detector, as well as providing 
the calibration information necessary to reconstruct events. 
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Slow-control data acquisition and execution of slow-control commands are 
carried out by the processor on the motherboard of the relevant electronics 
container (OM, LCM, SCM, or specialized instrumentation container). 

5.6 Calibration and positioning 

The pointing accuracy of the detector is determined largely by the overall timing 
accuracy of each event. This is a quadratic sum of terms due to 

• the precision with which the spatial positioning and orientation of the 
optical modules is known (o-geom); 

• the accuracy with which the arrival time of photons at the optical modules 
is measured (cjpmt); 

• the precision with which local timing of individual optical module signals 
can be synchronised with respect to each other {aaUgn) 

2 _ 2 . 2 I 2 

t pmt ' geom ' align 

Furthermore, in situ calibration of the optical module efficiency as a function 

of time is necessary in order to measure and correct changes in the response due 
to factors such as optical fouling. These issues are now discussed in detail. 

5.6.1 Positioning 

The reconstruction of the muon trajectory is based on the differences of the 
arrival times of the photons between optical modules. As such, it is sensitive 
to the distances between the optical modules. In order to avoid degrading the 
reconstruction, it is necessary to monitor the position of each optical module with 
a precision of 10 cm (light travels 22 cm per ns in water). The reconstruction of 
the muon trajectory and the determination of its energy also require knowledge 
of the optical module orientation with a precision of a few degrees. The precise 
absolute positioning of the whole detector has to be guaranteed in order to point 
to individual sources. 

To attain a suitable precision on the overall positioning accuracy constant 
monitoring of relative positions of the various detector elements with respect 
to absolutely positioned beacons is necessary. A full description of the relative 
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positioning equipment on the string and its performance during prototype string 
deployments has already been discussed in section |4.5.3 . 

The absolute positioning of the detector is performed by acoustic triangulation 
of low frequency acoustic beacons placed on the string bottom and a rangemeter 
on a surface ship, equipped with the Dynamic Global Positioning System (DGPS). 
A precision of ~ 1 m has already been demonstrated for the positioning of 
a prototype string, deployed at full immersion depth. Precise triangulation 
necessitates the knowledge of the sound velocity profile from the sea floor to the 
surface, which can be strongly distorted by thermal effects which depend on the 
season and current. This has to be measured by using a CTD profiler when the 
absolute positioning calibration is performed. 



5.6.2 Timing precision and calibration 

A master clock on shore, linked to Universal Time (UT) through the Global 
Positioning System network, permits to match events to transient astronomical 
phenomena such as gamma-ray bursts. The clock signal is distributed through the 
array network to each OM. 

A LED is located in each OM, facing upwards towards the PMTs further up 
in the string. It can be pulsed synchronously to the clock signal. This system will 
enable time synchronization between adjacent storeys in the same string. 

Calibration with external light sources will relay the initial calibration and 
monitor any possible drifts. The proposed system consists of optical beacons 
which contain high intensity pulsed light sources. They will illuminate several 
strings simultaneously. Ideally, the optical beacons should emit light at a 
wavelength as close as possible to that which is least attenuated in water (A ~ 
470 nm) and should distribute the light efficiently throughout the full solid angle 
with a pulse width of about 1 ns. Pulsed solid-state lasers and blue LEDs are the 



focus of detailed studies. Figure [5^ illustrates the range to be obtained from such 
an optical beacon at various wavelengths, as a function of pulse energy. 

Doubled Nd-YAG lasers which emit green light (A = 532 nm) can be passively 
Q-switched to give sub-nanosecond pulses. Lasers delivering pulse energies of 
~1 /iJ (~ 10^^ photons) and with a time spread (FWHM) of ~0.5 ns (see figure 
|S!9| , left) have been tested. These lasers are small, easy to operate and very 
robust, complying with the stringent deployment and installation requirements 
of the experiment. 

After a warming-up period of a few minutes Nd-YAG lasers show an output 



power stability better than 1% over several hours (see figure 5.9, right). 
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Figure 5.8: Distance in water reached by a pulse of light as a function of its energy for 
different wavelengths. At least 5 photoelectrons are required to be detected in the PMT. 
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Figure 5.9: Left: Time shape distribution (picoseconds) of a Nd-YAG laser pulse as 
measured by a streak camera. Right: Warming-up curve (minutes) of a Nd-YAG laser as 
seen by a photomultipher After a few minutes the output power stability is better than 
1% over several hours. 



Simulations have been performed to estimate the reach of optical beacons 
using lasers together with lambertian diffusers (~ cosO). For 1 /x J Nd-YAG lasers, 
distances in excess of 200 m can be reached. 



A LED pulser circuit based on that of Kapustinsky et al. [ |63| ] has been 
developed using recent GaAs LEDs from Nichia and Hewlett Packard which emit 
in the blue (A = 470 nm). The light pulse produced has a risetime of 2.0 ns and a 
duration of 4 ns FWHM as can be seen in figure |5.1(J| which was measured using 
the single photoelectron technique. 

The intensity of a single flash is 50 pJ (10^ photons). The timing jitter between 
the trigger pulse and the light flash is below 100 ps which indicates that, it is 
possible to enslave large numbers of such flasher modules together to produce 
bright flashes. 



5.6.3 Optical Module efficiency calibration 

The optical module efficiency may evolve with time due to modifications of 
the PMT characteristics (gain, quantum efficiency, etc.), or to a change in the 
transmission of light from the glass sphere surface to the PMT photocathode 
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Figure 5.10: LED time response 



(fouling, ageing of gel, etc.). As discussed previously, the K present in sea 
water, is a natural source of single photoelectron signals. The distribution of 
the collected charge from such events can therefore be used to monitor the gain 
of the photomultipliers. An in situ measurement of the OM efficiency using an 
external calibrated light source will be a convolution of three distinct effects: the 
fouling of the light source itself, variations of the light transmission in sea water, 
and the OM's intrinsic efficiency, i.e. the change in the response of the OM in 
units of single photo electrons. A relative calibration can, however, be done by 
comparing the response of several OMs to the high power light sources, which 
will be installed in the detector for the purpose of timing calibration. The light 
from these sources will arrive simultaneously at many OMs, allowing an accurate 
cross-check of the variation in amplitude response of the OMs. In the case of the 
Nd-YAG lasers, the accuracy of the measurement of the amplitude is limited by 
statistics, not by the source stability. Since the change in efficiency is expected 
to be a slow phenomenon, the statistical accumulation of calibration events will 
permit an accurate determination. In the same way, the LEDs on each storey 
pointing at the OMs further up the string can be used to measure the change in 
light transmission due to fouling and other phenomena. 



77 



To avoid problems caused by fouling of the light source, a ship-tethered 
calibration system is being considered for absolute efficiency calibration. This 
would consist of an isotropic calibrated light source enclosed in a Benthos 
sphere and a rangemeter to permit the accurate positioning of the source using 
the acoustic positioning system of the detector. An accuracy of a few tens of 
centimetres on distances of a few tens of metres between locations of the source 
in-between strings and the optical modules can be achieved and is satisfactory 
for our purpose. The optical module orientation with respect to the light source 
is determined from the tiltmeter and compass information. Sufficient fluxes of 
photons can be obtained by triggering simultaneously a large number of blue 
LEDs mounted inside a Benthos sphere. Absolute calibration of the optical 
modules requires that the attenuation length of light be known. This can be 
deduced directly from data taken with the light source at a few different positions 
with respect to a single optical module. 



5.7 Onshore data acquisition 



The object of the onshore data acquisition is to apply a level 3 filter to reduce 
this rate to a reasonable level for archiving on tape, maintain an experiment status 
database using the slow control information, and verify the integrity of the data. 

The first step of the on-shore processing is event building. Here, time-stamped 
data from various parts of the detector are assembled into events. Most of the 
triggers are caused by accidental coincidences, and a customised version of the 
off-line reconstruction will be used to filter these. 

The event building will associate both digital data from the optical modules 
and slow-control data in the same event. Fully built events will be used for 
feedback in controlling the detector, as well as for event display and data 
monitoring. A Unix-based event display has been developed and tested. The 
programme for monitoring of slow-control parameters and event histogramming 
has been written and tested for the prototype string connected to shore (see 
section P3[ ). It is based on the EPICS package and will be used in the 0.1 km^ 
detector. 
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Chapter 6 

Detector performance 



The neutrino sources discussed in the scientific programme fall into two distinct 
categories: 

• high energy neutrinos (above 1 TeV) from astrophysical sources; 

• low energy neutrinos (below a few hundred GeV) for oscillation studies and 
from the annihilation of neutralinos. 

The two classes behave differently in the detector — for example, muons from 
high energy neutrinos will pass right through the detector volume, whereas the low 
energy neutrinos may be partially contained — so both must be considered when 
optimising the detector design. Extensive simulation studies have been carried 
out on the detection of neutrinos from AGN, representing the first class of high 
energy astrophysical neutrinos, and on the measurement of neutrino oscillation 
parameters, which requires a good understanding of sub-TeV neutrinos. This 
chapter describes the event simulation and reconstruction, the optimisation and 
performance of the detector in both energy regimes, the results of the simulation 
studies, and a discussion of relevant systematic effects. 

6.1 Monte Carlo simulation tools 

This section describes the software tools used for event generation and detector 
simulation. Separate simulation packages are needed to treat signal events and 
various sources of background. 
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6.1.1 Event generation 



Neutrinos from astronomical sources have been simulated by generating the 
kinematic distributions of the neutrino interactions and taking into account the 
absorption and energy loss of the muons in the Earth. Events are generated with 
a flat spectrum in both the log of the neutrino energy and in the cosine of the 
angle of incidence, then weights are applied according to the flux of the sources. 
Events are retained for analysis if there is a muon track anywhere within a distance 
L = 140 m of the instrumented detector, as shown in figure ^7l\ 

The main background for upward-going astronomical neutrinos comes from 
neutrinos produced by cosmic rays interacting in the atmosphere. For atmospheric 
neutrinos interacting within the volume shown in figure |0| , tracks due to the 
hadronic showers, as well as those due to the muons, are generated using the 
LEPTO package [|4|]. For neutrinos interacting outside this volume, only the 
tracks due to the muons are simulated. 

Muons originating directly from the cosmic ray showers are simulated with 
a single muon component and a multi-muon component. High statistics samples 
of single muons have been generated with a parametrization of the flux at the 
detector depth [ |65| ] and also by propagating the sea level flux to the detector 
using the package PROPMU [^. The multi-muon component has been generated 
using the HEMAS package [ p7| ] with the muons propagated to the detector using 
PROPMU. At the present time, only relatively small samples of multi-muons have 
been produced. 



6.1.2 Detector simulation 

The response of the detector to the various types of physics events has been 
studied mainly using a version of the DADA program which originated in 
the BAIKAL collaboration. This program has been modified to improve the 
performance at high energies and to adapt it for the ANTARES detector. In 
DADA, the muon is tracked using GEANT, and light coming from the secondary 
particles generated is parameterized for energies above 0.5 GeV for electrons, 
above 1 GeV for photons, and above 10 GeV for hadrons. Less energetic 
secondary particles are taken into account by increasing the amount of Cherenkov 
light emitted by the muon (-1- 18% -i- 3.2% log (GeV)). A parametrization is 
used for the angular and longitudinal distribution of the Cherenkov light initiated 



by electromagnetic and hadronic cascades [68]. 



Checks on some of the simulation results have been made with a new program. 
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Figure 6.1: Scheme of the 'can' used for event generation. The distance L is taken as 
2.5 times the maximum absorption length for light. 



GEASIM, developed for ANTARES. Some of the simulations of multi-muons 
have also been performed with the KM3 program which, unlike the other 
programs, simulates the scattering of light in the water. 

If not specified otherwise, the simulation was performed with 10-inch 
photomultiplier tubes. All parameters used in the simulation are based on the 



measurements described in section 5.3.3, The angular response simulated is 



that shown in figure p3| , with a timing accuracy of 1.5 ns. For each optical 
module the time and amplitude of the pulse corresponding to an event is simulated 
taking into account the noise coming from the optical background, the gain of the 
photomultiplier and the effect of the electronics. All simulated events include 
random hits from '^°K in the sea water. 

Photons are scattered and absorbed by sea water. DADA includes the effect of 
absorption, with an absorption length of 55 m unless otherwise stated, but does not 
include scattering. For a distance of 45 m between source and optical modules, 
less than 10% of the photons are scattered, so the approximation of no scattering 
should have little effect on the signal. Nonetheless, the scattering must be included 
for the background studies. No additional smearing has been incorporated into the 
simulations for the relative positioning of the optical modules. 
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Figure 6.2: Average number of single photoelectron hits versus the number of 
coincidences for high energy muons, with and without optical background. The simulated 
background rate is 20 kHz. 

6.2 Track finding and reconstruction 

Before track reconstruction, a pattern recognition step is necessary to remove hits 
from the "^^K background. The hardware trigger will make requirements on the 
number of coincidences between two of the three optical modules in each storey 
of the detector. For the simulations of astrophysical performance a total of at least 
four coincidences on at least two different strings were required; the requirements 
for the oscillation study are discussed in section |^.4. 1| below. 



6.2.1 Pattern recognition 

Figure shows the number of hits with amplitude corresponding to one 
photoelectron, for a simulation of high energy muons with and without 
background. About 90% of single photoelectron hits come from . They cannot 
be removed by an amplitude cut without a large loss in signal hits. 

The method adopted to remove the ""^K hits is to start with a track fit which 



uses only coincidence hits. Figure |63| shows the angular error of this prefit when 
four coincidences are required. The prefit is used to define a road in time for the 
selection of single hits to be used in a further fit. The hits are kept if the time 
difference between measurement and expectation is between —150 and +50 ns 
and if the distance between the track and the optical module is smaller than 100 m. 
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Figure 6.3: Angular error for the prefit solution for high-energy events with at least 4 
coincidences on at least 2 strings, for an initial 1/ E'^ neutrino spectrum. 

This selection reduces the number of hits due to the optical background by a factor 
of 100, while keeping 46% of the single photoelectron hits associated with the 
muon tracks. 

6.2.2 Track reconstruction 

For each optical module, the arrival time of the Cherenkov light is 

{ti)o = to + {Li + (ii tan 6'c)/c 

where di, Li are defined in figure 6c is the Cherenkov angle, and to is the 
time at a reference point. Then (tj)o is smeared out by the time resolution 
(PMT transit time spread, positioning error, scattering) to give the pulse time 
ti. Photons from secondary particles will also be detected, with a distribution 
of arrival times that decreases exponentially with respect to (ti)o- Photons from 
the optical background arrive at random times. 

Figure ^3] shows the expected time distribution of the hits for different muon 
energies. This is the probability distribution function -P(At) which depends on the 
muon energy. In order to simplify the fitting procedure, the energy dependence 



83 



PMTi(Xi,yi,zi) 




Figure 6.4: Propagation of a muon and its associated Cherenkov light. 

is neglected. The probability distribution is evaluated either at the most-likely 
energy value or as an average over the energy range under study. 

First a preliminary estimation (prefit) of the track parameters is made, based 
on a linear fitting procedure [|70[]. Then the full fit is made, using a likelihood 
function (C) which is the product of the probability for each optical module i 
registered at time tj : 

i 

The maximisation of C by numerical methods yields the five track parameters 
X, y, z (or if:), 9, and 0, as well as a covariance matrix. This matrix provides the 
errors which will be used when applying selection criteria. 

6.3 Astrophysical neutrinos 

The principal interest of a deep-sea neutrino detector is the study of astronomical 
sources such as active galactic nuclei and gamma-ray bursters. These are 
expected to produce very high energy neutrinos yielding TeV-PeV muons which 
will not stop within the detector fiducial volume (uncontained events). This 
section presents the results of a simulation study for different models of neutrino 
production in AGN. 

6.3.1 Event selection 

After the track reconstruction, a selection is required in order to keep well- 
reconstructed events. The main selection criteria are: 
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Figure 6.5: Unnormalised probability density function P{/^t) for the arrival time of hits 
originating from muon tracks with four different energies: 2 TeV, 10 TeV, 50 TeV and 
250 TeV. 

• the errors on parameters of the reconstructed tracks have to be small: 

6^, 6y < 5m 6b <l° 5$ < 2° 6t <6 ns; 

• the distance between the initial muon track given by the prefit and the final 
track is required to be less than 150 m and the angle between them less than 
45°. 



These criteria have to be optimised for each physics channel. The effects 



of one set of criteria dedicated to neutrino astronomy are shown in figure |6^ . 
A dramatic reduction is observed in the tail of the distribution. These selection 
criteria will be used in the following to obtain the expected detector performances 
for diffuse sources of neutrinos. 



6.3.2 Optimisation for high energy events 

Given an overall conceptual design of a 0.1 km^ detector with about 1000 
optical modules, a number of detailed design parameters can be varied to 
optimise the performance of the detector. This process is not completely intuitive 
since the optimization procedure is dependent on the physical processes under 
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Figure 6.6: Angular error for the reconstruction witli and without selection (a) for the 
full range of errors, and (b) for the low error range only. 

consideration. The objective has been to define a detector with a good angular 
resolution, a large effective area and a good rejection against cosmic ray muons. It 
will be shown that the energy resolution is not heavily dependent on the geometry. 

The strings are arranged as shown in figure |6J| in order to avoid potential 
symmetries which could contribute to ambiguities in the reconstruction. 

The detector parameters that have been considered are: 

• the number of photomultipliers n in a local cluster (n = 2, 3, 4) and their 
orientation. 

n = 2 the PMTs look downwards (2) 

n = 3 the PMTs are at 45° below the horizon (3) 

n = 4 the PMTs are horizontal (4H) 

or the PMTs are at 45° below the horizon (4S) ; 

• the distance between two local clusters : dz = S, 12, 16, 20 m ; 

• the distance between two strings : dn = 60, 80, 100, 120 m. 

Other parameters such as PMT size and the optical properties of water affect 
the optimisation. It has been assumed that they are decoupled. The optimisation is 
made by changing them one at a time. Checks are made to ensure that the angular 
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Figure 6.7: Strings arrangement proposed to avoid symmetries. 



and the energy resolutions are not affected by the geometry, to first order. The 
uniformity of the detector response as a function of the zenith angle has also been 
considered. The following figures show the sensitivity of the results to the various 
detector geometry parameters. Optimisation for specific physics analyses will be 
considered later. 



The effect of varying the parameter dz is illustrated in figure The value of 
dz = 16 m has been chosen; the corresponding effective area is close to the value 
found for dz = Sm, and the zenith-angle response is more uniform. 

Figure ^|9| shows the effect of the number of photomultipliers in a storey 
(the total number of optical modules is kept constant). The geometry with three 
photomultipliers oriented at 45° below the horizontal is preferred to that with four 
horizontal photomultipliers on the same argument as above, namely uniformity in 
zenith angle. 

Figure |6.1U| shows the variation in effective area when modifying the 
parameter dn', the choice of dn depends on the energy range to be optimised. 
A 60 m spacing ensures the best effective area for energies up to some tens of 
TeV and becomes less efficient for higher energies. It should be noted that the 
effect of scattering in water has not been considered in this study. By taking this 
into account and thus having a smaller effective attenuation length it is expected 
that smaller string spacings will be favoured even at the higher energies. 

The configuration which was chosen for the following studies has 15 strings 
with an inter-string distance of 80 m. Each string is equipped with 21 storeys 
spaced vertically by 16 m, each storey consisting of 3 photomultipliers oriented at 
45° below the horizontal. 
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Figure 6.8: Effective detector area as a function of the muon energy when varying the 
vertical distance between two local clusters dz- The area is averaged over all muons 
coming from the lower hemisphere. 
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Figure 6.9: Effective detector area as a function of the muon energy when varying the 
number of photomultipliers (n) in a local cluster. The area is averaged over all muons 
coming from the lower hemisphere. 
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Figure 6.10: Effective area as function of the muon energy, varying the horizontal 
distance between two strings dn- The area is averaged over all muons coming from the 
lower hemisphere. 




Figure 6.11: Angular resolution for a neutrino spectrum, for at = 1.3 ns. 
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Figure 6.12: Median angles as a function of log{Ey). The total angle, the physical angle 
and the reconstruction error are indicated. 



6.3.3 Angular resolution 

The angular resolution depends on the algorithms used, on the selection applied 
to the reconstructed events and on the timing accuracy which was discussed in 



section EJi 



The angular resolution is shown in figure |6.1 1| for at = 1.3 ns. Half of the 
events from a simulated point source with a l/E"^ spectrum fall inside a circular 
region of radius 0.2°. One hemisphere can be effectively divided into 200,000 
pixels. 

It is worth noting that above 10 TeV the angular resolution is not dominated 
by the physical angle between the neutrino and the real muon, but by the 



reconstruction errors (see figure 5.12) 



6.3.4 Energy resolution 

For energies above 1 TeV, the muon energy loss is dominated by catastrophic 
energy loss from bremsstrahlung and pair production, both of which increase 
with energy. The amount of light produced by the track is proportional to the 
energy loss and this allows an estimation of the muon energy. An estimator x 
which is roughly proportional to the energy is formed using the number of hits. 
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their observed pulse-height amplitudes Oj, and the amplitudes 0^(0) expected for 
a minimum-ionizing particle : 



X = Nhit I - 1 ) 

Since the energy information comes from stochastic processes, a truncated mean 
is used to decrease the fluctuations. Hits are selected for the energy estimation 
depending on whether they are lower-energy hits (smaller pulse heights) or higher- 
energy hits (larger pulse heights). Only hits which are within 4 ns of the arrival 
time expected from the reconstructed track (|tj — (ti)o| < 4 ns) are candidates 
for use in the estimators. Hits are included in the calculation of the low- and 
high-energy estimators if they satisfy one of the following criteria : 

• 0.1 < < 100 for the low-energy estimator, 

• 10 < < 1000 for the high-energy estimator. 



The number of hits satisfying the selection criteria is shown in figure |6.13K a) for 



a lower-energy muon (2 TeV) and in figure |6.13K b) for a higher-energy muon 



(50 TeV). The low-energy estimator is used if more than 10% of the hits are 
accepted for that estimation; otherwise the high-energy estimator is used. 

The average value of log(x) as a function of \og(E^) is shown in figure |6.14 



for events selected for the low- and high-energy estimations. The linear portions 
of these curves are used to estimate the muon energy between 1 TeV and 10 PeV. 
The following parametrizations are obtained : 



1.56 + 1.04 X log(x) for the low-energy estimator 
1.26 + 1.02 X log(a;) for the high-energy estimator 



Figure |6.15| shows the reconstructed energy spectrum for simulated muons 



with energies from 1 TeV to 1 PeV. The resolution on logE is 0.40 at 1 TeV, 
decreasing to 0.33 at 100 TeV; this implies a muon energy resolution AE/E of 
about a factor 3 below 10 TeV, improving to a factor 2 above 10 TeV. 

In figure |6.16|, it can be seen that taking into account the energy resolution 



enables a correct energy spectrum to be unfolded. By placing a cut on the 
reconstructed energy (figure |6.17[ ) it is expected that the contamination of 
atmospheric neutrinos can be controlled. 
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Figure 6.13: Number of hits used for the energy evaluation (a) at 2 TeV, using the low 
energy muon estimator, and (b) at 50 TeV, using the high energy estimator. 
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Figure 6.14: Parametrization of the low energy (a) and high energy (b) estimators. 
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Figure 6.15: The reconstructed muon energy for different true energies. 




Figure 6.16: True and reconstructed spectra for (a) muons coming from atmospheric 
neutrinos and (b) muons coming from AGN (NMB model). 
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Figure 6.17: Reconstructed spectra of muons coming from (a) generic AGN and 
(b) blazars, compared to muons coming from atmospheric neutrinos (ATM). 

6.3.5 The diffuse neutrino flux from active galactic nuclei 

The integrated neutrino flux from a large population of distant AGN generates a 
diffuse neutrino flux (similar to the diffuse X-ray background) even if individual 
sources are not detectable. Calculations of this diffuse flux have been done for 
four models of neutrino production by AGN: 

• SDSS [[7T|]: generic AGN model in which the interactions of the protons 
with the material of the accretion disk dominate; 

• NMB [[7^: generic AGN model in which the dominant interaction is with 
the ambient radiation field; 

• PRO [|73|]: blazar model where the protons are accelerated in the innermost 
region of the jets and interact principally with accretion disk radiation; 

• MRLA-B [[7^: lower and upper limits, respectively, of a blazar model in 
which the protons are accelerated in the jets and interact primarily with 
synchrotron radiation produced by relativistic electrons in the jets. 

Other calculations of diffuse neutrino sources are described in reference [Q]. 
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Figure |6.18| shows the effective area for the events which trigger the detector, 
for those which are reconstructed and those which satisfy the selection criteria 
(see section |OTTt ). 
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Figure 6.18: Effective area for triggering, reconstructed and selected events. The areas 
are averaged for muons coming from the lower hemisphere. 



The corresponding counting rates of accepted events are summarised in 
table ^]T] for different diffuse sources (atmospheric neutrinos and AGNs) in the 
optimised detector with different (true) muon energy thresholds. For the same 
sources, the expected number of events above the reconstructed muon energy 
threshold is given in table |^ . 

The sensitivity is model dependent, and is affected by the rejection of cosmic 
ray muons and by the measurement of the muon energy. Beyond 10 TeV, if no 
event rate above that for atmospheric neutrinos is observed (this corresponds to ~ 
40 events per year), then a minimum detectable l/E'^ flux at the 5 a level of 



dK 



> 6 - 10 



4[ \ 
UTeV; 



is obtained (assuming that the down-going cosmic muons can be successfully 
rejected). 
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Model 


Minimum E*™^ 


lOGeV 


1 TeV 


10 TeV 


100 TeV 


Atmospheric 
ATM 


2725±146 


544±44 


18±2 


0.34±0.05 


Generic AGN model 

SDSS 

NMB 


354±13 
591±15 


334±12 
454±13 


260±12 
206±9 


132±10 
53±4 


Blazars 
PRO [73] 
MRLA 
MRLB 0] 


45±2 
34±1 
37±1 


43±2 
18.1±0.7 
21.1±0.8 


35±2 
4.0±0.2 
6.6±0.4 


21±2 
0.52±0.04 
2.4±0.3 



Table 6.1: Diffuse neutrino flux from AGN. The table shows the number of events per 
year with true muon energy greater than E^^'^'^ which are reconstructed as upward-going 
in the optimised detector The errors are statistical. 



Model 


Minimum E'^,'^" 


10 TeV 


100 TeV 


Atmospheric 
ATM 1^] 


68±13 


0.8±0.1 


Generic AGN models 
SDSS [71] 
NMB [[72] 


251±12 
217±9 


134±10 

64±4 


Blazars 
PRO [0] 
MRLB [g 


34±2 
7.8±0.4 


21±2 
2.6±0.3 



Table 6.2: Diffuse neutrino flux from AGN: the number of events per year with 
reconstructed muon energy greater than E^^^ which are accepted as upward-going in the 
optimised detector. The errors are statistical. 
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6.3.6 Detecting astrophysical point sources 



Due to the extremely low angle between the muon and the parent neutrino and 
to the good quality of the muon direction measurement, the atmospheric neutrino 
background contaminating each individual source can be reduced to a very low 
level by selecting very small angular regions of the sky. For a point-like source, 
visible 50% of the time, evidence for observation will be possible with 6 events 
above 1 TeV (3 in the central pixel). This corresponds to a fake probability smaller 
than 7%. 

The sensitivity for an individual source would be: 



-2 



-1 



The corresponding luminosity in protons emitted at any energy can be derived 
with the usual assumption that 5% of the proton energy is converted into neutrino 
energy. This gives a sensitivity of about 10^° W (about 1/10 the Eddington limit 
for a solar mass) for an object located at a distance of lOkpc. 

Another way to estimate the ANTARES potential for the detection of point- 
like sources is to calculate the number of events expected for a given flux. A flux 
estimate can be obtained by using the measured low-energy gamma-ray fluxes, 
assuming i) that the low energy gamma-rays are of hadronic origin and ii) that 
the emitted gamma-rays have a differential energy spectrum With these 

assumptions the muon neutrino flux is about 40% of the flux of gammas at the 



production source. Using the second EGRET catalog QTSQ for sources measured 
during the PI 2 period, the derived neutrino flux has been extrapolated to the 
energies where the neutrino detector is sensitive. 

The detection of the most luminous individual sources could require more 
than one year of data acquisition. Nevertheless, a statistically significant effect 
could be detected in one year by adding the contributions of all the extra-galactic 
sources. The expected number of events is between ten and one hundred per year 
(depending on the value of the differential spectral index used: 2 or 2.2) to be 
compared to a total background of about three events per year. 

The sensitivity to gamma ray bursts is improved by the detection of events 
correlated in space and in time with optical information. In this case, the 
background is negligible. With the flux given in reference [ [T0| ] and the variations 
of the effective detector area shown in figure |6.18| , the rate expected is ~ 10 signal 



events per year with a background rate of less than 0.001 events per year. 
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6.3.7 Background rejection 



One critical item is the rejection of down-going muons from the normal cosmic 
ray flux. At a depth of 2300 m in water, down-going muons correspond to a 
flux 1.5x10^ times higher than that for the atmospheric neutrinos. Therefore, a 
rejection factor of more than 1.5x 10^ is required to have a signal-to-background 
ratio exceeding 1 . This is less critical for point-like sources where the rejection 
required is decreased by the number of pixels defined in the sky (about 200,000). 

The rejection is, in principle, easy to study but in practice the results are 
limited by the number of triggering cosmic ray events that can be generated. 
Simulated single muons and multi-muons reaching the detector have been studied. 



as discussed in section 6.1. The simulation statistics are summarised in table 5.3 



Approximately three events per second will satisfy the combined trigger and 
filter conditions, of which 38% are single muons and the remaining 62% are multi- 
muons (half of them coming from primaries of more than 200 TeV, very few from 
initial energies smaller than 2 TeV). 

For the potentially most dangerous events (multi-muons coming from a very 
high energy primary cosmic), the equivalent of 4 days of data has been generated. 
Clearly more statistics are required for an effective background estimate; however, 
the existing sample took 80 CPU-days of a single processor at the Lyon computer 
centre. A factor of 3 improvement in simulation speed is expected as a result 
of using the KM3 program instead of DADA for the detector simulation, and a 
further factor of 2 gain is possible if the fiducial volume is reduced. 

From all available simulations, no cosmic ray muon or multi-muon is accepted 
as an upward-going muon after reconstruction and quality tests (see figure |6.19D . 
Since the most dangerous event class will be the multi-muon events initiated at 
primary energies > 200 TeV, this implies, with the current statistics, a rejection 
factor Nfake/Natm.u < 1/3 at 90% confidence level. AMANDA and BAIKAL do 
not experience problems due to this background. 



6.3.8 Systematics on detector performance 

The performance of the detector depends on the values of various simulation 
parameters relating to the properties of the photomultipliers, the transparency of 
the water and the number of strings. In the actual detector these values are likely to 
be different from the preliminary numbers used in the simulations, so it is essential 
to study the effect of such variations on the performance of the detector. 

As an example of the effects to be expected, the effective area of the detector 
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mono-muons 




multi-muons : 








2-20 TeV 


20-200 TeV 


> 200 TeV 


Flux at 2300 m 
depth (cm~^s~^) 


10-^ 


0.4-10"^ 


1 • io-« 


0.4 • 10"^ 




85% 


15% 






22% 


56% 


22% 


No. of simulated 

primary 
interactions 




141 • 10^ 


28 • 10^ 


12 • 10^ 


No. of simulated 
events reaching 
the detector 

level 


0.3 • 10^ 


22 000 


0.7 • 10^ 


6- 10^ 


No. of 

triggering 

events 


27800 ± 167 


230±70 


15 093 ± 183 


384 746 ± 620 


corresponding 
real time 


< 1 TeV 1.5 days 
1-10 TeV 0.9 days 
> 10 TeV 20 days 


30 minutes 


4.5 hours 


3.9 days 


no of triggering 
events per hour 


4698±2U 


460±im 


5290±41 


4074±7 




38% 


62% 






6% 


42% 


52% 














12 522 ±262 
100% 



Table 6.3: Summary of the simulation statistics. For the multi-muons, is the primary 
energy. 
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Figure 6.19: Number of downward-going muons as a function of the reconstructed zenith 
angle. 



as a function of the muon energy, previously shown in figure |6.18| , has been 
recalculated using the same detector layout with parameters based on more recent 
work: the most recent optical module angular dependence measurement, an 
absorption length of 47 rather than 55 m, an angular resolution of 0.3° and an 
overall timing error of 2.5 ns. The result is shown in figure p.20 . 

A more quantitative result is obtained by studying each parameter separately, 
beginning with photomultiplier parameters. The final values of these are not yet 
known since the choice of phototubes has not yet been made. Ten-inch tubes are 
currently used in the simulations since they are functional and readily available. 
However, photomultipliers with a larger photocathode are expected in the near 
future. 

Going from a 10-inch to a 12-inch PMT, the potential rate for al/E'^ spectrum 
increases by 40%, and at increases by 15%. Therefore, the overall gain on S/ \fB 
would be 10%. 



Figure |6.21| shows how an increase of the absolute efficiency of 20% and a 
change in the angular response of the optical module affect the effective detector 
area. These parameters need to be measured accurately for a good understanding 
of the systematic errors. 

The choice of phototube also affects the timing error ou and consequently the 
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Figure 6.20: Effective area for the optimised detector with the measured parameters. 



angular resolution ag: ag varies roughly as as shown in figure |5.22| . The 
8-inch tubes used in OM tests to date have apmt = 11 ns, which combines with 
an alignment error of about 0.5 ns to give an overall timing error of 1.3 ns and a 
corresponding angular resolution of 0.2°. For the 10-inch PMT, the contribution 
of the PMT is larger (apmt = 1-5 ns); if, in addition, the alignment error were as 
bad as 2 ns, the angular resolution would be 0.3°. 

For a given string spacing, the effective area of the detector clearly depends on 
the number of strings deployed. Figure |6.23| shows the relative effective area of 8 
and 45 string detectors, normalised to the standard 15-string detector, as a function 
of E^. The effective area grows linearly with the number of optical modules at 
low energy, but the effect decreases above 10 TeV. 

The angular resolution improves when the number of strings increases: for 
eight strings ag is 20% greater than the 15-string value, and for 45 strings it is 15% 
less. Better angular resolution leads to improved sensitivity to point-like sources. 
In practice, a larger detector would not simply be an extension of the present array: 
operating experience with the 0. 1 km^ detector and additional optimisation studies 
would be required to develop an appropriate design. 

The effective area is also a function of the water transparency. In the sea water 
at the detector site, effective attenuation lengths of 40-60 m have been measured. 



as discussed in section 4.2.3. The effect of varying the absorption length in the 
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Figure 6.21: Effective area sensitivity to PMT parameters: (a) standard parameters; (b) 
increase quantum efficiency by 20%; (c) change of the angular response. 
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Figure 6.22: The angular resolution as a function of the overall time resolution. 
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Figure 6.23: Relative effective area versus the logarithm of the muon energy for 45 
strings (triangles) and 8 strings (crosses), normalised to 15 strings. 
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Figure 6.24: Effective areas for absorption lengths of 40 and 55 m. 



103 



simulation studies from 40 m to 55 m is shown in figure p.24\ It should be noted 
that the reconstruction and selection criteria were optimised for 55 m, and may 
not be optimal for 40 m. 

The scattering effect gives fewer delayed photons (5% at 25 m, 10% at 45 m) 
than those produced in catastrophic muon energy loss processes, so the angular 
resolution should not be affected. 

6.4 Neutrino oscillations 

Neutrino interactions below a few hundred GeV can be studied in ANTARES 
because the neutrino energy can be estimated from the range of contained muons 
observed in the detector volume. These events can be used to study additional 
physics channels such as neutrino oscillations and neutralino annihilations. This 
section describes the simulations that have been carried out to understand the 
sensitivity to oscillations of atmospheric neutrinos. 

The principal source of neutrinos for neutrino oscillation studies is the decay 
of charged pions produced by cosmic rays interacting in the Earth's atmosphere. 
Charged-current interactions of producing upward-going muons are selected 
so that the background from muons produced in the atmosphere is absorbed by 
the Earth. 

The simulation programs used for the astronomical neutrino studies are also 
used for the oscillation studies. The reconstruction algorithms, on the other hand, 
are somewhat different. This chapter will describe the current status relative to 
the simulation and reconstruction of the signal and background for the neutrino 
oscillations studies, and the determination of the oscillation parameters. 

6.4.1 Track reconstruction 

For the study of neutrino oscillations, events are reconstructed in two classes 
depending on the number of strings with coincidences. If at least two strings have 
coincidences, the event is classified as 'multi-string', while if only one string has 
coincidences, the event is classified as 'single-string'. However, only 7% of the 
multi-string events have only two coincidences and none of the well reconstructed 
single-string events has less than two coincidences. 

For the multi-string neutrino oscillation events, a prefit similar to the one 
described earlier is used, but with somewhat tighter cuts. For the single-string 
events, a method is used which estimates the angle of the track to the string using 
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Figure 6.25: Illustration of the event containment estimator a. 



npe- 



the time differences between adjacent storeys; then three hypotheses are tried for 
the distance of the track from the string and the one with the best fit likelihood is 
used as the start for the final fit. 

For the neutrino oscillation studies, the muon range is used to measure the 
muon energy from which the neutrino energy is estimated. Events must be 
selected which are at least partially contained within the detector volume. For 
the single-string events the containment is ensured by requiring the reconstructed 
muon track to have a zenith angle within 20° of the vertical and requiring the 
two highest and two lowest storeys of the string not to have fired in the event. A 
partial-containment condition is used for the multi-string events, where tracks are 
used with zenith angles as large as 80°. Figure |05| illustrates the algorithm used. 
Once the track is reconstructed, the points, where the first (A) and the last (B) 
detected Cherenkov photons were associated to the track, are located. Then the 
track is extrapolated backwards and a test is made to see if any photomultipliers 
could have seen that part of the track before the point A. If none could, then the 
track is clearly not contained. 

The estimator of the track containment is the sum of the amplitudes of the hits 
that should have been seen on these optical modules but were not seen. A simple 
model is used to calculate these amplitudes. A minimum-ionizing particle seen 
at normal incidence by an optical module at a distance of 1 m yields an average 
of 55 photoelectrons [|6T|]. The sum of the amplitudes for a minimum ionizing 
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Figure 6.26: Distribution of a„pe for contained and non-contained events. 



106 



track emitting Cherenkov light at an angle Oj with respect to optical module i at 
a minimum distance di from the track (figure |6.25[ ) is: 



Here /(A0) accounts for the angular response of the optical module. Figure [06 
shows the distributions of this quantity for contained and non-contained events; a 
cut on anpe at a value of 50 provides an efficient separation between the two sets 
of events. 



6.4.2 Optimisation for oscillations 

The observation of atmospheric neutrinos for energies below a few hundred 
GeV can be used to study neutrino oscillations for values of the neutrino mass 
difference Am^ between 10^'^ and 10^^ eV^. The detector configuration proposed 
above for neutrino astronomy has been evaluated for neutrino oscillations. Other 
studies have been carried out to see if the optimal geometry could be different 
for neutrino oscillations and astronomy. Simulations have been made with five 
different geometries, varying the number of optical modules per storey and 
their orientation (pointing horizontally, down at 45°, or vertically down). The 
conclusion is that three optical modules per storey, pointing down at 45°, is also 
optimal for neutrino oscillations. A denser vertical spacing is optimal for the 
oscillation studies; therefore, the geometry proposed here has four strings with 
8 m vertical spacing and nine strings with 16 m spacing. The studies performed 
so far are not exhaustive, and a better compromise may be found using the same 
total number of strings and optical modules. A horizontal spacing between strings 
of 60 m is optimal both for astronomy below 100 TeV and for oscillations. 



6.4.3 Event selection 

After reconstruction of an event, the quality of the reconstruction is estimated 
using the fitted likelihood and the error estimates provided by the likelihood fit. 
Another variable used to discriminate against badly fitted tracks is the angular 
distance between the track found in the prefit and the track found in the final fit. 
All tracks are required to have at least 10 hit optical modules. The mean number of 
hits associated with the tracks is 30 for multi-string tracks and 16 for single-string 
tracks. The mean number of hit coincidences (two out of three hits in a storey) 
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A© (deg) 



Figure 6.27: Angular resolution. The shaded histogram is the difference in angle 
between the direction and the true fi direction and the open histogram is the complete 
reconstruction error between the f direction and the reconstructed fi direction. 



is 7.2 for multi-string events and 5.3 for single-string events. Selected events 
are required to satisfy the containment and zenith angle conditions described in 
section |6.4.1| . 

After reconstruction and containment cuts, the number of charged-current z/^ 
events expected, in the absence of neutrino oscillations, is 2400 for three years of 
data taking with the proposed detector. About 90% of these events are partially- 
contained multi-string events; the remaining 10% are fully-contained single-string 
events. 



6.4.4 Angular and energy resolution 



The neutrino angular resolution for the selected events is shown in figure |6.27 



The experimental angular resolution is small compared to the intrinsic difference 
between the neutrino and muon directions. The median angular resolution from 
measurement errors is 0.4° for the multi-string events and 1.8° for the single-string 
events. 

The correlation between the reconstructed muon energy and its true energy 
is presented in figure |6.28| for the two types of events. The single-string events 
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Measured Muon Energy vs. True Muon Energy (GeV) 
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Figure 6.28: Correlation between reconstructed and true muon energy for tracks 
reconstructed as multi-string and single-string. 

are mainly concentrated at low energy. The true muon energy of the multi- 
string events extends above 100 GeV because the muon end-point is not required 
to be inside the detector volume. The intrinsic muon energy resolution of the 
ANTARES detector is adequate for the oscillation measurements. Nonetheless, 
attempts will be made to improve the resolution on the neutrino energy by 
estimating the energy of the individual hadronic showers. 

The spatial resolution for the starting point (vertex) and stopping point for 
the multi-string events is shown in figure |6.29| . The multi-string events can be 



reconstructed in all their spatial dimensions with reasonable accuracy, but only the 
vertical coordinate is well measured for the single-string events. The distribution 
of the reconstructed zenith angle and the energy for the muons from the partially- 
contained events described above are shown in figures |6.30| (a) and (b). 



6.4.5 Backgrounds 

The reconstructed zenith angle and the muon energy for events from 
interactions which pass all the cuts as partially-contained events are shown in 



figure |6.30| (e) and (f). The events amount to 4% of the signal events; their 
angular distribution is similar to the signal but they are at lower mean energies. 
Figure |6.30| (c) and (d) show the distributions for the z/^ events which are truly 
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Figure 6.29: Resolution on the track starting point and stopping point for multi-string 
events. 
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from interactions outside the detector, but which pass the cuts as contained events. 
These events are 20% of the signal events. Their reconstructed energy has a 
mean of 55 GeV, but their real energy is higher - a mean of 140 GeV - since they 
come from interactions beyond the detector which pass the containment cuts due 
to inefficiencies in the light detection. 

Another potential source of background is the upward-going muons from the 
decay of pions and kaons from neutral-current interactions of neutrinos of all 
flavours. This background is small, however, less than 4% of the charged- 
current interactions. Moreover, since these events have lower energies due to the 
missing final-state neutrinos, the acceptance will be low, and the number of fake 
events coming from this background can be neglected. 

Background from muons formed in the atmosphere are more time consuming 
to generate, due to the immense rate, than the background from neutrinos. The rate 
of muons arriving at the detector at angles truly below the horizon is completely 
dominated by i/^ interactions (figure |3TT1). But high statistics samples are needed to 
measure the rate of badly reconstructed down-going muons. With the simulations 
performed to date, no events from atmospheric showers pass the selection cuts, 
but the statistics simulated correspond to only two days of data. This sets a limit 
on the possible number of atmospheric muons mis-identified as z/^t interactions, 
but higher statistics are needed for a full evaluation. 



6.4.6 Sensitivity to oscillations 

An analysis has been carried out using simulated interactions to explore the 
sensitivity of the proposed experiment to the neutrino oscillation parameters Am^ 
and sin^26'. The analysis is based on partially-contained muons, and the muon 
energy is taken as the visible energy in the detector. The number of signal 
events was calculated using the presently known values for the flux of atmospheric 
neutrinos and the expected experimental acceptance. The neutrino flight distance 
Li, is related to the neutrino zenith angle 9z by L^^ ~ Lq cosOz, where Lq is the 
diameter of the Earth (taken as a perfect sphere of diameter 12740 km). 

The histogram in figure |6.31| shows the distribution of accepted events as a 
function of E^/(L^/ Lq), where is the energy and is the neutrino flight 
distance calculated using the zenith angle of the reconstructed muon track, in 
place of the zenith angle of the neutrino. The histogram represents a simulation 
of three years of operation without neutrino oscillations. The 'data' points in 



figure 5.31 represent the same simulated events, but with each event weighted by 
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Figure 6.30: Distributions of reconstructed muon angle and energy: (a),(b) partially- 
contained signal Vn events; (c),(d) non-contained events and (e),(f) background 
events. 
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Figure 6.31: Simulated number of events in 3 years with no oscillations (histogram) and 
with oscillations (points) for sin^ 26 = 1.0 and Am^ = 0.0035 eV^. 



the probability: 

P = 1 - sm^2e sm^{1.27Am^L^/E^), 

using the true neutrino values of and L^. 

To explore the Am^ and sin^2^^ parameter space, a series of three-year 
simulations were carried out, each one producing a distribution of 'data' points 
like the points shown in figure |6.31| . A large number of such simulations were 



carried out for each point in parameter space. Each of these distributions was 
compared to a high-statistics simulation without oscillations. An example of the 
ratio oscillations/no-oscillations for one of these three-year simulations is shown 
in figure |6.32. 



These simulations have been used to find the regions of parameter space that 
could be excluded if no oscillations are found in the data. To remove the effect 
of the absolute normalization, the no-oscillations distribution is multiplied by a 
factor which is varied as a free parameter for each of the three-year simulations. 
Then the for the hypothesis that all of the 'data' points in the spectrum are 
compatible with no oscillations is calculated. If the probability is greater than 
10%, then the experiment is said to be compatible with no neutrino oscillations 
for the given parameters. Conversely, in the region of parameter space where 
the probability is less than 10%, evidence for oscillations could be claimed with 
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Figure 6.32: Simulated spectrum of a typical experiment for 3 years of data taking 

showing the ratio of data corresponding to oscillations with sm^26 = 1.0 and Am? = 
0.0035 eV^, to Monte Carlo expectation with no oscillations. The horizontal line is the 
average level used to test the statistical significance of the oscillations and the solid line is 
the fit used to extract the oscillation parameters. 
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Figure 6.33: For three years of data taking, neutrino oscillations could be excluded at 
90% confidence level in 80% of simulated experiments in the region to the right of the 
curve. 



90% confidence. Due to statistical fluctuations there is a spread in probabilities 
between different experiments for the same parameters. The exclusion plot shown 
in figure |6.33| gives the region where 80% of experiments could give an exclusion 



at the 90% confidence level. 

The simulations are also used to estimate the accuracy with which the 
oscillation parameters can be extracted from the data. Distributions such as those 
shown in figure |6.32| are fitted with a damped oscillation function: 



Pfit = 1 - sin' 29[1 /2 - (1/2 - sin' y)e~^^^y^' 



where y = (1.27Am'L^/ E^) / a/1 + P and / is a damping parameter. In the 
analysis this damping parameter is obtained from a fit to a high-statistics Monte 
Carlo spectrum for each set of oscillation parameters. Then the 'data' distributions 
from the three-year simulations are fitted with / fixed and Am' and sin'2^^ 
free. Figure |6.34| displays the fitted oscillation parameters compared to the true 



generated parameters. A good correlation between the fitted parameters and the 
true values is seen over a large range of the parameters but there is a bias in the 
fitted value of Am' because the muon energy is smaller than the neutrino energy. 
In the experiment, this bias would be corrected using the Monte Carlo. 
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Figure 6.34: Results of the fits to neutrino oscillation parameters: (a) fitted and true 
values of Am^ for fixed sin^ 26 (true) = 1, and (b) fitted and true values of sin^20 for 
fixed Am^ (true) = 0.0035 eV^. 



Figure |6.35| shows the region of parameter space where both Am^ and siv?26 
can be measured with an error less than 33% in more than 80% of experiments. 
It can be seen that the experiment can make accurate measurements for 0.002 < 
Am^ < 0.008 and sin^26' > 0.8. The measurement region is shifted up in Am^ 
compared to the exclusion region shown in figure |6.33| because of the low energy 
cut-off in the acceptance. A range ofE/L values around the minimum of the 
survival probability is needed to find the precise position of the minimum and 
make an accurate determination of Am^, whereas a measurement of the number 
of events at the calculated position of the minimum is sufficient for exclusion. 
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Figure 6.35: For three years of data taking, the oscillation parameters Am^, sin^20 can 
be measured with fractional error less than 33% in more than 80% of experiments in the 
region to the right of the curve. 

6.4.7 Systematics 

The present studies have shown that the proposed ANTARES detector can make 
useful measurements with the statistics expected in three years. It is important to 
note, however, that the present analysis is exploratory and can be refined in many 
ways. The quasi-elastic events (estimated to be about 15% of the charged- 
current events) were not included in these simulations. These events have a 
better correlation between the muon energy and the incident neutrino energy than 
the deep-inelastic events used here. Inclusion of these events will improve the 
experimental sensitivity to low Am^. On the other hand, the sensitivity will suffer 
in the real experiment due to systematics resulting from variations in the detector 
acceptance and assumptions about the neutrino flux and its energy and angular 
dependence and about backgrounds. 

The event reconstruction used at present has separate procedures for multi- 
string and single-string events. In the future, a uniform event reconstruction will 
be developed to avoid the systematic errors caused at the edge of the acceptances. 

The backgrounds have so far not been taken into account in the analysis. 
However, making a known background subtraction with the magnitude of 
backgrounds presently estimated would have little impact on the results. The 
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backgrounds from z/g and from non-contained z/^ are related to the atmospheric 
flux and are well simulated, so a background subtraction can be made. The present 
simulations indicate the backgrounds from muons originating in the atmosphere 
are small. These simulations must be refined by including light scattering in 
the water, and extended to higher statistics. The detector will have a significant 
acceptance for events above the horizontal, and these events can be used to verify 
the muon background simulations for events below the horizontal. 

In the analysis presented above, it was assumed that the shape of the neutrino 
flux was known as a function of energy and zenith angle, but that the overall 
rate was unknown. A different analysis can be performed which uses the non- 
contained events to normalise the contained events. This method would 
have similar statistical errors to the method presented here, but very different 
systematics. 

The knowledge of the flux of atmospheric neutrinos in the energy range from 
10 GeV to 1 TeV is of crucial importance for the study of neutrino oscillations, 
because observed deviations from the expected flux could be interpreted as a 
possible oscillation signal. Several calculations of this flux exist [ fT^ , [77], [7S], [79| ]. 
The comparison of the total flux shows variations of the order of 20% between the 
different calculations. The energy dependence of the flux is also slightly different 
in the various approaches. The difficulties in obtaining consistent results come 
from assumptions about the primary cosmic ray flux as well as from the lack of 
knowledge concerning the development of atmospheric hadron showers. 

A compilation of present experimental data on upward-going muons is shown 
in figure |6.36| , taken from Fogli et al. Only the zenith- angle dependence is 
available. The experimental data show large fluctuations, but they are generally 
compatible. Nonetheless, even the combined data shown in figure |6.36| (b) have 
at least a 20% bin-by-bin uncertainty. This is likely to be the limiting systematic 
on neutrino oscillation measurements by ANTARES, and studies of theoretical 
calculations and archive data are under way to reduce it as much as possible. 
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Comporison of upword — going muon data 
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Figure 6.36: Comparison of all upward-going muon spectra, rescaled to a common 
threshold of 1 GeV, for (a) data from individual experiments, and (b) combined data. 
The figures are taken from Fogli et al. [8C]. 
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Chapter 7 

Complementary techniques 



A number of related projects for the construction of neutrino telescopes in either 
ice or water are currently underway. The status of these projects is reviewed in this 
chapter. The case of AMANDA is given particular consideration, as the project is 
presently making rapid progress. 



7.1 Deep ice detector: AMANDA 



AMANDA [|8T|] was installed in several stages in the deep Antarctic ice at the 
South Pole. The detector has a rather simple design: the analogue signals 
of optical modules (OMs) are directly transmitted to surface through 2 km 
cables. In 1993-94, 80 OMs of AMANDA-A were deployed on four strings 
at depths of 810 to 1000 m; however, measurements of the ice transparency at 
those depths showed that the light scattering was unacceptable for maintaining 
good directional information. Subsequently, in the Austral summer 1995-96, 
86 OMs on four strings (AMANDA-B4), were deployed at depths of 1500- 
2000 m. The AMANDA-B4 array was mainly used for calibration, and testing 
track reconstruction and cuts on the real data. The scattering length for Cherenkov 
light in AMANDA-B depths was found to be two orders of magnitude larger 
than in shallower depths. In the following year, six new strings have completed 
AMANDA-B, which now consists of 302 PMTs on 10 strings. First analysis 
of the data from 85 days of effective lifetime resulted 16 up-going neutrino 
candidates, where 21 up-going atmospheric neutrinos were expected from Monte 
Carlo simulations. 

During the 1997-98 season, the construction of AMANDA-II, which will have 
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an effective area of 30,000 has begun. Three 1 .2 km long strings were installed 
with modules that range from 1150 to 2350 metres in depth. In the 1999/2000 
season 6 additional strings will be deployed to complete the AMANDA-II array. 
The new strings will test a variety of new techniques, such as the usage of 10" 
photomultiplier tubes, and 25-50 digital optical modules. 

AMANDA is also proposing a km scale array (ICECUBE) of 5000 OMs on 
80 strings to be deployed by 2007-08 at 1.4 to 2.4 km depths in the Antarctic ice. 

7.1.1 AMANDA compared to ANTARES 

Apart from their technological differences, the main difference between 
ANTARES and AMANDA is the use of water versus ice as the detector medium. 
Ice has the advantages of great stability (no currents, hence no movement of 
deployed strings), no living organisms to cause optical fouling or bioluminescence 
and low noise environment due to absence of ^°K. The advantages of water are a 
much longer scattering length, as measured by ANTARES test deployments, and 
flexibility in experimental design (strings can be retrieved and re-deployed). 

The difference in scattering length between ice and water is substantial: 
AMANDA measures 24 ± 2 m for ice at 1700 m whereas ANTARES test 
deployments indicate a scattering length in excess of 200 m. This leads to a 
significant improvement in the angular resolution of ANTARES. An accuracy of 
0.3° is achievable. 

7.1.2 Astrophysical sources 

From the point of view of sky coverage, a detector at the mid-northern latitudes is 
an ideal complement to the South Pole location of AMANDA. Neutrino telescopes 
look downwards, and are therefore sensitive to sources which are not circumpolar 
at their location: AMANDA 'sees' the northern celestial hemisphere (sources 
with positive declination 6). ANTARES, if located at latitude 42° 50' N, wiU 
see all sources with declination 5 < +47° 10': when combined with AMANDA 
this provides complete sky coverage, with overlap for confirmation and cross- 
calibration in the region 0° < 5 < 47°. It is worth noting that the Galactic centre, 
a region of interest for neutrino astronomy because of its high level of activity, 
is in the southern hemisphere and is therefore visible to ANTARES but not to 
AMANDA. 

The low noise environment of AMANDA makes it more suitable for 
the detection of low-energy (MeV) neutrinos such as those associated with 
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supemovae. However, neither AMANDA nor ANTARES is really designed 
to compete with dedicated low-energy neutrino detectors such as Super- 
Kamiokande SNO |Q or the proposed OMNIS for supernova 

observations. 



7.1.3 Neutrino oscillations and neutralino searches 

Measurement of neutrino oscillations by the method described above (effectively, 
measuring the energy spectrum) requires a low energy threshold, closely spaced 
optical modules and good background rejection. As detailed in chapter |^, 
ANTARES should be capable of exploring the region of the Super-Kamiokande 



signal [|28[] in a few years of running. It may be more difficult to do the energy 
measurement by range in a more diffusive medium, such as the South Pole ice. 

As discussed in chapter 2, neutralino searches by neutrino telescopes rely on 
the detection of neutrinos produced when neutralinos gravitationally captured in 
the centre of the Earth or the Sun annihilate. The angular distribution of muons 
produced by neutrinos from neutralino annihilation in the Earth's core can be used 
to infer the neutralino mass [Q, therefore, a good angular resolution is needed 
to make best use of any signal. The calculated width of the neutrino-induced 
muon distribution ranges from about 2.5° for high to ~ 7.5° for ~ 80 
GeV/c^, so an angular resolution of better than 1° is highly desirable. A pair of 
neutralinos does not annihilate directly into a pair of neutrinos, so the neutrino 
energy spectrum peaks at around one-third to one-half of the neutralino mass. 
Therefore, a fairly low energy threshold is useful. 

In ANTARES both the angular resolution and the low energy threshold can 
be achieved by using the analysis technique proposed for the neutrino oscillation 
search. The AMANDA search [ pS] ] uses a conventional multi-string analysis and 
loses sensitivity below a neutralino mass of about 100 GeV/c^. 



7.1.4 Complementarity 

Detection of a new phenomenon by a single experiment is never entirely 
satisfactory. The fact that ANTARES and AMANDA are broadly comparable 
experiments with similar capabilities in many areas promises rapid, independent 
confirmation of any new discovery, avoiding the long period of uncertainty 
associated with, for example, the solar neutrino deficit. Moreover, as shown 
above, ANTARES and AMANDA are also in many respects complementary: 
AMANDA has lower background noise and is already operational, ANTARES has 
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a longer scattering length, hence better angular resolution, and can be optimised 
for specific physics targets. The two together provide excellent all-sky coverage 
for astrophysical point sources. Overall, the combination of ANTARES and 
AMANDA addresses all the physics goals of high-energy neutrino telescopes 
better than either one could do alone. 

7.2 Deep water detectors 

7.2.1 Baikal 

The Baikal [ |8^ Neutrino Telescope is situated in the Siberian Lake Baikal, at a 
depth of 1.1 km. It consists of eight strings supporting 192 OMs. This project 
started in April 1993 by deploying 36 OMs in three half-strings. The system was 
recovered after 300 days of running time. The main problem was the biofouling 
deposited on the glass spheres, mainly on the upwards facing optical modules. 
In April 1996, four full strings carrying 96 OMs were deployed (NT-96). Only 
two layers of OMs were facing up. In winter 1997, 144 OMs were deployed, and 
the detector was completed in 1998 with 192 OMs. In the analysis [ [S7| ] of 70 
days of data taken with NT-96 array 2.0 x 10'' events were reconstructed, nine of 
which were up-going. The expected number of atmospheric neutrinos from MC 
studies was 8.7, and the angular distribution of the events was consistent with that 
of atmospheric neutrinos. The effective area of the detector had a maximum value 
of 400 m^ for this analysis. A different method was developed specifically to look 
for vertically up-going events, for which an effective area of about 1000 m^ was 
reached. This analysis found 4 events with zenith angles > 150° where 3.7 events 
were expected. From this result Baikal set a 90% C.L. upper limit of 1.1 x 10~^^ 
cm^^s^^ on the muon flux from the centre of the Earth for energies > 10 GeV. 

7.2.2 NESTOR 

NESTOR is planned to be installed in a 3 800 m deep sea site off Pylos 
(Greece). The project consists of a tower of 12 hexagonal floors of 16 m in radius 
supporting 168 OMs in total. The deployment of aluminium and of titanium floors 
has been tested at 2 500 m depth. Digital transmission of the signals is being 
developed. 
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7.3 Long baseline experiments 



The anomalous results seen in various atmospheric neutrino experiments could, 
in principle, be explained either by z/^ — oscillations or — oscillations, 
even though the latest Super- Kamiokande results favour z/^ — z/^. The — z/g 
explanation has essentially been ruled out by the results of the CHOOZ [j8^] 
experiment. The CHOOZ collaboration observed electron neutrinos at a distance 
of one km from a reactor complex by observing inverse /3-decay in a large 
heavily- shielded gadolinium-loaded liquid scintillator target. The disappearance 
of electron neutrinos as they oscillated to muon neutrinos would have resulted 
in a distortion of the observed electron spectrum from inverse beta decay. No 
distortion was seen, and the combination of the long baseline and low background 
meant that the entire parameter space for — z/g oscillations allowed by the 
existing atmospheric neutrino experiments was ruled out. A second long-baseline 
reactor experiment is underway at Palo Verde in the United States to verify this 
result. 

Explanations for the observed atmospheric neutrino anomaly based on z/^ — 
Vt oscillations could be tested by means of accelerator-based long-baseline 
experiments. In these experiments, an almost pure muon neutrino beam is 
produced from the decay in flight of a pion beam at an accelerator facility and 
directed towards a distant underground laboratory. These muon neutrinos produce 
muons in the detector by charged-current (CC) interactions. The first signature of 
neutrino oscillations would be a change in the total muon rate or a distortion in 
the observed energy spectrum caused by the energy-dependent disappearance of 
muon neutrinos. Of course any neutrino, with or without oscillations, can produce 
neutral-current (NC) interactions in the detector, so a second signature for neutrino 
oscillations would be a measured NC/CC ratio different than that expected for a 
pure muon neutrino beam. A third signature would be the appearance of r leptons 
coming from CC interactions of r neutrinos resulting from oscillations. There are 
currently three projects: K2K 0, MINOS [0, and CERN/LNGS. The potential 
advantage of ANTARES is in the large fiducial mass of the detector which will 
allow the measurement of atmospheric neutrinos above 5 GeV. 



7.3.1 K2K 

The 12 GeV primary proton beam at KEK is used to produce a muon neutrino 
beam directed at the Super-Kamiokande detector 250 km away. It first passes 
through two near detectors on the KEK site. As the energy of the resulting 
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neutrinos is too low to produce r leptons in the final state, the experiment searches 
for muon neutrino disappearance. 

The K2K experiment sent its first neutrinos towards the Super- Kamiokande 
detector in March 1999. The sensitivity of the experiment is limited by the 
luminosity of the beam. However, the experiment has the advantage of a 
considerable head start over the other projects. Figure pTT] displays the sensitivity 
after three years of running. 



7.3.2 MINOS 

The MINOS experiment will consist of a specially built beam-line directing 
neutrinos from the 120 GeV proton beam of the new Main Injector at 
Fermi National Accelerator Laboratory (FNAL) to a 5.4 kton magnetised iron 
calorimeter to be built in the Soudan mine, 730 km away in northern Minnesota. 
The existing Soudan II detector will form an additional far detector with a smaller 
fiducial mass but finer-grained calorimetry. There will be a near detector on the 
FNAL site. 

The initial version of the detector is limited to muon neutrino disappearance 
signatures. The group plans to upgrade their detector in the future, with a hybrid 
emulsion detector which would be capable of imaging r decays. The predicted 



sensitivity of the experiment is shown on figure 7.1 



7.3.3 CERN/LNGS 

There is much interest in constructing a third long-baseline facility using a beam 
from CERN directed towards the Gran Sasso Laboratory. This idea is proceeding 
to the point of formal proposals to be submitted to the respective laboratories in 
autumn 1999. Because they emphasise r appearance signatures, the CERN/LNGS 
projects may be best suited for precise measurements. 



7.4 Gamma ray telescopes and Air Shower arrays 

Gamma ray telescopes and neutrino telescopes are obviously complementary, 
though they may probe the same sources. Purely electromagnetic sources would 
not be seen in neutrinos, but hadronic sources will give both neutrinos and 
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photons. Neutrinos are less likely to be absorbed at the production site or during 
the transport (photons above 10-100 TeV are absorbed by 77 interactions). 

The Pierre Auger project is designed to measure very high energy cosmic 
rays (charged and neutral). Specifically, the experiment will also detect the 
showers caused by ultra-high-energy neutrinos (10^^ eV) travelling horizontally. 
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Figure 7.1: Expected sensitivity for K2K and MINOS compared with 
Kamiokande and Super- Kamiokande signals. 
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Chapter 8 

Cost and schedule 



The cost of a detector consisting of 4 strings of 123 optical modules and 9 
strings of 63 optical modules, for an overall number of 13 strings and 1059 optical 



modules, has been estimated. Table |0] gives a summary of the cost of these two 
types of strings. 



Cost of one string 


63 OM 


123 OM 




kFF 


kFF 


Acoustics: 
4 rangemeters, 
concentrator, interface 


500 


500 


Optical modules 


1200 


2340 


Top of string 


120 


120 


Bottom of the string 


390 


390 


Qualification, 
transport, assembling 


35 


50 


elementary segments 


825 


1610 


Electronics cost, 
bottom of the string 


150 


175 


Container 


1050 


2050 


Total 


4270 


7235 



Table 8.1: Cost estimates for the two types of strings. 
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Funding is already available for the construction of the first string: the 
'demonstrator". The costs for the further 12 strings are listed below: 

• Strings: 

4 with 123 OM, 8 with 63 OM 70 MFF 
including 10% contingency 

• Electro-optical cables: 

2 cables / 1 2 fibres / 40 km 9 MFF 

• Deployment operations: 

80 days costed, 10-15 needed at sea 10 MFF 

• Sea floor interconnections and junction box: 4 MFF 

• Main power supply and transformer 2 MFF 

• Electronic cost for trigger, DAQ on the shore: 2 MFF 

• Total cost: 97 MFF 



These estimates are based on current experience. Large scale production could 
lead to lower costs and enable extra strings to be built. 

Construction of a 1000 optical module detector would proceed in two stages. 
The first stage would be a ~ 500 OM detector, costing 48 MFF, to be built as soon 
as possible and enabling physics to start. For this stage the electronics would be 
the critical path item. 

The proposed schedule would be to build: 

1. Seven strings (3 high density, 4 low density) in 1999-2000-2001, to be 
deployed in 2002; 

2. The remaining strings to be ready to deploy in 2003. 
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Chapter 9 
Conclusion 



In the summer of 1996, the ANTARES collaboration started to study the 
feasibility of mounting a neutrino telescope in the Mediterranean Sea near 
Toulon. The first phase of the project was to measure parameters of the marine 
environment and to acquire and develop the technology necessary to perform an 
undersea experiment. 

This phase is now completed. A location has been chosen and explored using 
the submarine Nautile. The topography has been mapped and the sea current 
measured during periods covering two years. The rate and nature of sedimentation 
at the site has been measured in different seasons, and the loss of transparency of 
the optical modules was found to be less than 2% per year. The transmission 
properties at the site have been measured for blue light. The absorption length 
has been determined to be of the order of 50 m, whereas the effect of scattering is 
found to be negligible. Studies of the optical background in the sea water due to 
natural radioactive decays and biological organisms indicate that these rates can 
be dealt with. This background will lead to an inefficiency of the detector of less 
than 5%. 

The technology required for the experiment has been attained with the 
indispensable help of IFREMER and INSU, as well as through developments 
with industry. The prototype string has proven our ability to deploy a detector 
string and the submarine cable connections have been demonstrated. The design 
of mechanics and electronics for a full detector is well advanced and most essential 
elements have been illustrated. 

This proposal has presented the extensive scientific programme possible with 
an undersea neutrino telescope. The same apparatus could detect cosmic neutrinos 
in the range from TeV to PeV with the ability to measure the diffuse neutrino 
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flux from active galactic nuclei, make searches for dark matter in the form of 
relic neutralinos, and measure neutrino oscillation parameters. Evidence for point 
sources of neutrinos could be obtained by integrating the neutrino flux from 
the positions of known candidate sources such as AGN and gamma-ray bursts. 
Simulations have been performed to show the physics potential of the detector 
and optimize its geometry. 

The next step is the proposed construction of a detector consisting of 13 strings 
with about 1000 photomultipliers in an array covering an area of 0.1 km^. The 
detector will be completely modular and can be deployed in the sea in stages. The 
present schedule calls for 7 strings to be deployed in 2002 with continual tests 
and partial deployments at intermediate stages. Physics measurements could start 
with this initial detector. The full 13 string detector should be operational in 2003. 
The cost of the proposed detector is estimated at 15 MEuros (97 MFF). 

The ANTARES project has now grown to a size of over 100 members 
from five European countries. While still actively searching and attracting new 
collaborators, the present institutes have the manpower and resources required to 
build and operate a detector of the size and complexity proposed. 

The collaboration between particle physicists and astrophysicists has a long 
and fruitful history. It has given rise to entirely new ways of looking at the 
Universe, such as X-ray and gamma-ray satellite-based telescopes, and recently 
ground-based gamma ray telescopes extending the field of astronomy above 
10 TeV. There is every reason to believe that neutrino astronomy will open a new 
window on the Universe at least as exciting and intriguing. 

Some of the goals outlined in the scientific programme cannot be achieved 
on a realistic time scale with a 0.1 km^ detector. For example, to observe 
most astrophysical point sources a detector on the square kilometre scale will 
be required. The experience gained in constructing and operating the detector 
described here will be essential for the success of a km-scale project. 
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